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ABSTRACT 


This  report  presents  the  results  of  the  past  year's  work  in  a  continuing 
program  to  investigate  basic  processes  in  thermionic  energy  conversion 
important  to  a  thermionic  nuclear  power  plant  for  naval  applications.  The 
previous  work  was  reported  in  AI-6?99»  "First  Summary  Report  of  Basic 
Research  in  Thermionic  Conversion  Processes."  The  subjects  discussed  in 
the  present  report  are: 

Statistical  Mechanics  of  Cesium  Adsorption 

Space  Charge  Analysis  for  Low  Pressure  Thermionic 
Diodes 

Emission  Requirements  for  Removal  of  Space  Charge 
Barriers 

Unignited  Mode  of  Thermionic  Converters 

Interpretation  of  Volt- Ampere  Characteristics 

Vaporization  and  Deposition  at  Cesium  Covered 
Surfaces 

Cesium  Purification 
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I.  INTRODUCTION 

This  research  program  is  directed  toward  development  Of  an  understand¬ 
ing  of  the  fundamental  processes  occurring  in  thermionic  converters  and 
correlation  of  these  processes  to  develop  an  understanding  of  overall 
converter  performance 4 

The  previous  phases  of  the  program  were  mostly  concerned  with  emission 
processes;  these  results  were  reported  in  AI-6799,  "First  summary  fieport  of 
Basic  Research  in  Thermionic  Conversion  Processes, "  issued  in  November  1961, 
In  that  report  the  ideal  converter  performance  was  determined  and  the 
details  of  the  fundamental  emission  processes  were  analyzed  and  related  to 
the  requirements  for  ideal  converter  operation.  This  thereby  provided  an 
extensive  analysis  of  emission-limited  converter  operation. 

The  continuation  of  the  program ,  the  results  of  which  are  included  in 
the  present  report,  extends  some  of  the  theories  of  the  emission  processes, 
but  emphasizes  the  processes  occurring  in  the  interelectrode  space  and  the 
relation  and  importance  of  the  emission  characteristics  on  the  transport 
effects.  The  work  of  this  year,  therefore,  represents  a  logical  extension 
of  the  earlier  work  in  that  the  knowledge  of  the  emission  processes  is 
coupled  with  studies  of  the  transport  processes  in  order  to  effect  a  more 
complete  understanding  of  thermionic  energy  conversion.  The  results  are 
summarized  in  the  following  section.  The  remainder  of  the  report  consists 
of  detailed  articles  on  problems  investigated  during  the  past  year. 

Predominantly,  the  conventions  used  in  this  report  are  those  proposed 
for  the  1962  Colorado  Springs  Symposium  on  Thermionic  Power  Conversion, 

When  the  discussion  is  directly  related  to  gaseous  discharges  or  vacuum 
electronics,  however,  and  when  extensive  reference  is  made  to  published 
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works  in  those  fields,  the  terminology  and  conventions  of  those  fields  are 
retained  to  enable  the  reader  to  relate  conveniently  the  present  and 
previous  work.  For  example,  in  sections  IV  and  V  the  power  quadrant  is 
plotted  to  the  left  while  in  Sections  VI  and  VII  it  is  plotted  to  the  right . 
also ,  the  terms  cathode  and  anode  are  used  instead  of  emitter  and  collector 
when  discussing  previous  literature  which  use  the  former  terms.  The  terms 
emitter  and  collector  are  used  whenever  thermionic  converters  are  discussed.. 
When  subscripts  to  symbols  are  to  indicate  an  electrode,  upper  case  letters 
are  used;  when  subscripts  or  superscripts  are  to  indicate  a  particle  species 
lower  case  letters  are  used. 
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II.  SUMMARY  OF  RESULTS 

Some  of  the  work  of  the  past  year  has  been  a  continuation  of  the  study 
of  emission  phenomena  and  the  emission-limited  thermionic  converter.  The 
results  of  this  study  are  discussed  in  Section  III.  Results  of  the  studies 
of  the  transport  processes  and  the  interaction  of  emission  and  transport 
effects  are  discussed  in  Sections  IV,  V,  VI  and  VII.  Two  other  experiments 
were  performed  during  this  year.  One  of  the  experiments  involved  vapori¬ 
zation  and  deposition  at  cesium-covered  surfaces;  the  other  concerned  chem¬ 
ical  analyses  of  commercially  available  cesium.  The  results  of  these  two 
experiments  are  reported  in  Sections  VIII  and  IX.  Summaries  of  these 
Sections  are  as  follows : 

STATISTICAL  MECHANICS  OF  CESIUM  ADSORPTION 

A  statistical  mechanical  formalism  is  developed  for  the  adsorption  of 
cesium  on  metals.  The  adsorbed  cesium  is  assumed  to  exist  in  both  ionic  and 
atomic  states.  This  formalism  describes  more  accurately  the  approach  used 
in  the  first  ONR  report ,  Al-6799*  The  earlier  treatment  assumed  that 
adsorption  took  place  only  as  ions.  The  present  treatment  yields  general 
expressions  for  the  atom,  ion  and  electron  evaporation  rates  (v  ,  v.  and  v  ) 
with  both  aaions  and  adatoms  present.  Previously ,  v  ,  but  not  v  or  V. , 
was  derived  under  the  same  assumptions;  only  adions  were  assumed  and  the 
derivation  of  and  lay  partly  outside  the  statistical  mechanical 
formalism. 

The  inclusion  of  the  effects  of  adatoms  has  shown  that  a  major  effect 
in  explaining  the  maximum  in  the  work  function  change  may  be  due  to  the 
polari zability  of  the  adatoms.  Furthermore,  at  large  coverages  where  this 
theory  is  inadequate ,  the  lateral  bonding  of  the  cesium  and  possible  accom- 
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panying  effects  such  as  changes  in  polarizability  may  also  be  extremely 
important. 

SPACE  CHARGE  ANALYSIS  FOR  LOW  PRESSURE  THERMIONIC  DIODES 

A  space  charge  analysis  of  a  collisionless ,  plane  electrode  thermionic 
converter  has  been  performed.  With  both  ions  and  electrons  emitted,  the 
potential  distributions  can  assume  a  variety  of  forms.  This  report  completely 
systematizes  these  solutions.  The  space  charge  problem  is  first  considered 
for  infinite  spacing,  and  the  potential  distributions  are  characterized  by 
the  collector  potential  and  the  potential  at  which  the  electric  field  first 
vanishes*  For  a  given  emitted  ion-to-eiectron  ratio,  the  diode  ion  and 
electron  currents  are  easily  extracted  as  a  function  of  collector  potential, 
and  thus  the  volt-ampere  curve  can  be  plotted.  The  resulting  volt-ampere 
characteristics  extend  continuously  from  minus  infinity  to  plus  infinity, 
so  that  the  load  line  always  yields  an  operating  point.  The  present  analysis 
concludes  with  a  discussion  of  the  solutions  for  finite  spacing.  They  are 
compared  with  the  infinite  Spacing  solutions  and  numerical  results  are  given* 
An  example  of  bi-stability  is  shown,  and  the  origin  of  oscillations  briefly 
considered* 

EMISSION  REQUIREMENTS  FOR  REMOVAL  OF  SPACE  CHARGE  BARRIERS 

Potential  distributions  in  low-pressure  converters  have  been  analyzed 
theoretically  to  determine  under  what  conditions  they  will  be  operating  in 
retarding-field,  space-charge-limited  or  saturated-emission  modes.  The 
necessary  calculations  have  been  made  so  that  the  particular  mode  can  be 
determined  from  a  graph  when  the  voltage,  emitter  temperature,  spacing  and 
emission  8  of  the  converter  are  given*  The  same  data  allow  one  to  determine 
the  i  required  for  saturated  diode  current  for  given  diode  voltage,  emitter 
temperature  and  spacing.  One  simple  result  of  the  calculations  is  that  f,or 


AI-7979 

4 


@>1.27  a  space  charge  limitation  cannot  exist  in  a  diode  for  any  current, 
spacing  or  voltage.  Another  result  is  that  for  small  applied  voltages 
(eV>2kt)  an  emission  3»0.5  insures  that  no  space  Charge  limitations  will 
occur  for  any  value  of  electron  emission  current, 

UNIGNITED  MODE  OF  THERMIONIC  CONVERTERS 

The  volt-ampere  characteristics  of  thermionic  converters  in  the  unig¬ 
nited  mode  are  sensitive  to  changes  in  emitter  temperature,  cesium  pressure 
and  spacing.  The  dependence  of  the  characteristics  on  these  parameters  has 
been  measured  so  that  some  of  the  basic  processes  occurring  in  converters 
could  be  studied  under  conditions  not  complicated  by  the  additional  mechan¬ 
isms  in  the  ignited  mode.  Several  features  were  observed  electrically  and 
visually  which  related  to  previous  work  with  externally  heated  hot  cathode 
discharges.  Examples  are  the  "anode  glow"  and  "ball  of  fire"  modes  studied 
by  Malter,  Johnson  and  Webster.  The  first  ignited  mode  studied  by  Daly  and 
Emelens,  and  Pengelly  and  Wright  was  also  observed  in  thermionic  converters. 
Fateh  effects  were  also  noted  as  evidenced  by  variations  in  the  glow  in  the 
interelectrode  spacing  and  in  the  dependence  of  the  low  pressure  character¬ 
istics  on  voltage  and  spacing.  This  work  provides  the  experimental  basis 
for  the  theoretical  analyses  of  Section  V. 

INTERPRETATION  OF  VOLT-AMPERE  CHARACTERISTICS 

The  set  of  volt-ampere  curves  in  Section  VI  have  been  interpreted 
in  terms  of  surface  ionization  modes,  sheath  ionization  modes  and  volume 
ionization  modes.  The  emphasis  in  the  present  paper  has  been  on  the  first 
two  modes.  In  the  unignited,  surface  ionization  mode  two  limiting  modes 
are  studied,  a  space  charge  mode  and  a  plasma  mode.  The  space  charge  mode 
is  treated  by  space  charge  theory  and  the  plasma  mode  by  transport  theory. 
The  space  charge  mode  occurs  predominantly  for  spacing  which  is  small  com- 


AI-7979 

5 


pared  to  the  mean  free  path  and  vice-versa  for  the  plasma  mode.  The  close 
spacing  volt-ampere  curves  were  interpreted  by  the  space  Charge  analysis 
developed  by  Langmuir.  A  patch  effect  appears  to  be  indicated.  A  transport 
theory  is  developed  for  the  plasma  mode  and  the  importance  of  diffusion  as 
a  current-limiting  mechanism  is  experimentally  observed  in  the  dependence  of 
the  saturation  current  on  spacing.  This  observation  means  that  theories  of 
the  unignited  mode  which  assume  a  uniform  plasma  and  thus  no  density  gradient 
are  incorrect.  The  development  of  a  collector  sheath  sufficiently  large  to 
provide  sheath  ionization  is  observed  in  the  departure  from  the  apparent 
saturation  current  to  increased  currents  as  more  ions  are  provided.  A 
theory  with  sheath  ionization  is  Outlined  but  not  carried  out.  One  very 
important  feature  of  the  sheath  ionization  mode  is  its  initiation  at  sheath 
voltages  appreciably  less  than  the  cesium  ionization  potential.  Further 
confirmation  of  the  correctness  of  this  conclusion  would  indicate  a  two- 
stage  process  as  has  been  suggested  for  the  ignited  mode. 

VAPOHX ZAT ION  AND  DEPOSITION  AT  CESIUM  COVERED  SURFACES 

The  effect  of  adsorbed  impurities  (particularly  cesium)  on  vaporization 
rates  and  condensation  coefficients  has  been  discussed,  at  cesium  coverages , 
cesium  pressures ,  and  electrode  temperatures  of  interest  in  thermionic 
energy  conversion,  it  was  found  that  vaporized  emitter  material  has  only  a 
0.0005  to  0.001 3  probability  for  sticking  to  the  collector.  No  effect  with 
cesium  coverages  was  observed  and  close  correspondence  was  found  for  a  case 
where  no  back  scatter  occurred  to  a  case  where  normal  converter  back  scatter 
occurred.  It  was  concluded,  therefore ,  that  the  low  condensation  coeffi¬ 
cient  is  the  dominating  factor  determining  material  transport  in  converters. 
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CESIUM  PURIFICATION 

Since  it  has  been  shown  experimentally  that  impurities  or  additives  to 
cesium  vapor  have  a  profound  effect  on  converter  performance,  a  program  was 
initiated  to  examine  the  purity  of  the  commercially  available  cesium.  At 
the  beginning  of  the  program  considerable  uncertainty  existed  regarding  both 
the  purity  of  the  metal  and  the  accuracy  of  the  chemical  analyses  of  the 
metal,  analytical  techniques  were  established  for  examining  cesium  samples 
from  several  suppliers.  It  was  found  that  the  quality  of  commercially 
available  cesium  has  continually  improved  and  that  in  most  cases  the 
accompanying  analyses  are  adequate. 
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III.  STATISTICAL  MECHANICS  OF  CESIUM  ADSORPTION 

C.  Warner,  III 


INTRODUCTION 

A  statistical  mechanical  formalism  is  developed  in  this  study  for  the 

adsorption  of  cesium  on  metals.  The  adsorbed  cesium  is  assumed  to  exist  in 

both  ionic  and  atomic  states.  This  formalism  describes  more  accurately  the 

approach  used  in  the  earlier  paper  "Correlation  of  Electron,  Ion  and  Atom 

Emission  Energies."^"  Another  paper  in  that  report  "Statistical  Mechanical 

2 

Treatment  of  Surface  Ionization"  assumed  that  adsorption  took  place  only 
as  ions.  The  present  treatment  yields  general  expressions  for  the  atom, 
ion  and  electron  evaporation  rates  (v  ,  v.  and  v  )  with  both  adions  and 

3-1  6 

adatoms  present.  Previously,^"  v.  was  derived  but  not  v  or  v.  under  the 

©  &  1 

same  assumptions.  In  Reference  2,  only  adions  were  assumed  and  the  deri¬ 
vation  of  v,  and  v-  lay  partly  outside  the  statistical  mechanical  formalism, 
a  x 

In  Reference  3,  an  approximation  to  the  energy  of  ionic  evaporation  was 

made  with  adions  only.  The  approximation  considered  the  adions  always  to 

be  in  their  final  configuration  as  the  test  adion  was  being  adsorbed. 

The  results  of  this  study  are  compared  with  those  of  Taylor  and 

4  12  5 

Langmuir,  and  the  previous  ONR  work.  ’  ’  Finally,  inadequacies  in  the 
analysis  are  discussed.  In  particular,  the  difficulties  of  calculation  at 
high  coverages  are  discussed.  A  high  coverage  calculation  which  includes 
depolarization  of  adatoms  is  shown  to  alter  appreciably  the  work  function 
predictions. 

DEVELOPMENT  OF  GENERAL  EXPRESSIONS  FOR  DEPARTURE  RATES 

Consider  the  closed  system  at  constant  temperature  T  and  pressure  P 
shown  in  Figure  1,  and  let  there  be  N  cesium  atoms ,  N.  cesium  ions  and  N 

9  A 
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Figure  1.  Cesium- [Metal]  Systems 
at  Constant  T  and  P 


electrons  in  the  gas  phase  in  equilibrium  with  N  electrons  in  the  metal 

em 

and  Nga  cesium  atoms  and  Ng^  cesium  ions  on  its  surface*  For  phase  equili¬ 
brium  ,  thermodynamic  laws  require  that  the  corresponding  chemical  potentials 


must  satisfy  these  equations: 

•  V  *a  =  V’  wi  * 

and  for  chemical  equilibrium,  the  following  equation 
%  "*i  +  %  * 

The  usual  expressions  for  ugm,  and  are 


...  (1,2,3) 


...  (4) 


tt  m  kT  In  _e 

2V  mkT 


-  •* 0  [-J 


...  (5) 


...  (6) 


U,  *  kT  in  fi  2jA^ 
V  MkT 


u  =  kT  In  _a  /  2t£; 

2V  MkT 


...  (7) 


...  (8) 
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for  volume  V,  Ionization  potential  X*  electron  mass  m,  cesium  mass  M, 
Planck  constant  h  =  2tt~,  and  Boltzmann  constant  k.  The  work  function  0  in 
Equation  6  is  assumed  to  depend  only  on  the  ionic  coverage,  Ng^/o/  and  is 
given  by  electrostatic  theory. 


0  =  0  -  00  =  0  -  2TT  Nsi  M 

O  O  —7— 


...  (9) 


where  M  is  the  electric  dipole  moment  of  the  adion  and  its  image  and  J  is 
the  surface  area,  - — - - - - 


Equations  1,  5  and  6  imply 


fel  /kT  _  4TTmk3  T3  exp  (  - 
V  V  2irm  “h3 


) 


...  (10) 


which  is  just  Richardson's  equation  if  the  electron  departure  rate  is 
equated  to  the  electron  arrival  (which  is  the  left  hand  side  of  Equation  10 
by  kinetic  theory). 

Equations  4,  5»  7  and  8  imply  the  Saha  equation 


Na  ^e  /  2tt_T  ha  exp  (e  I/kT) 
N±  '  V  1  mkT| 


The  atom  and  ion  departure  rates,  and  v^,  follow  from  Equations  2 
and  3  just  as  the  Richardson  equation  follows  from  (1).  It  is  necessary  in 


obtaining  the  Richardson  equation  from  (10/  and  the  departure  rates  va  and 
v±  derived  below  in  Equations  24  and  25,  that  the  incident  particles  are 
not  reflected  but  all  adsorbed;  otherwise  the  departure  rate  must  be 
equated  to  that  fraction  of  the  incident  particles  which  are  adsorbed. 
This  assumption  implies  that  the  adsorption  process  does  not  have  any 
activation  energy. 

The  chemical  potentials  u  and  u  .  are  not  given  in  the  literature 
and  must  be  calculated  by  statistical  mechanics: 
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ZERO  OF  ENERGY 


Figure  2.  Energy  Curves  for  Cesium  Atoms  and  Ions 
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S3 


M 

1 

sg*a 

?/?//// 


METAL 


Figure  3»  Breakup  of  Energy  of  Ionic  Adsorption 

Corresponding  to  Three  Steps,  S, ,  S_,  S, 

1  <  j 


Figure  4,  Potential  Distribution  for  Finite 
Assembly  of  Dipoles 
(Points  X^  and  Xg  Correspond  to  Points 
Shorn  in  Figure  3) 
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X* 


M>  =  -  kT  S  In  Q 

88  3F“ 

sa 


U  .  =  -  kT  _A  -  In  Q 
81  W~ 

si 


T,/t  N 


si 


T,*/,  N 


sa 


where  Q  is  the  classical  partition  function  for  the  particles  in  the  surface 
phase  p 


N 


sN  +N  . 
sa  sa. 


<k  = 


3N  3N 


N  1  N  ,1  h 
sa  si 


sa 


si 


d3x,d3p.e 


i=l 


and  H  is  the  complete  Hamiltonian  of  the  adsorbed  particles.  The  following 
expression  is  considered  to  approximate  the  partition  function  Q: 


N  l  N  ,  J 
sa  si 


sa 


N 


sa 


N  +N  .  N  . 

sa  si  ^  si 


N„.<*„..+eIVkT 
.  sa  ao 
e 


eA  <t> 


/kT 


,N  .] 

sa  si 


(14) 


where  the  factor  2  is  a  result  of  the  spin  degeneracy  of  the  adatom;  qfi  = 
kT/hf  and  q.  =  kT/hf .  with  f  and  f .  being  the  frequencies  associated  with 
the  simple  harmonic  approximation  for  the  atoms  and  ions  in  the  direction 
normal  to  the  metal  surface;  Ash/  Y§nMkT  and  arises  from  the  momentum 
integrations  in  the  plane  of  the  surface ;  the  exponential  factors  are 
associated  with  the  zero  of  energy  which  has  been  taken  to  be  zero  for  ions 
or  electrons  at  rest  in  the  gas  phase.  In  Figure  2  the  surface  atomic 
state  differs  from  the  zero  of  energy  by  0  +  el;  the  ionic  energy  shift 

SO 

can  be  understood  from  Figures  3  and  4.  Figure  3  shows  the  container 
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(Figure  1)  on  top  of  a  metal  of  infinite  extent*  The  ion  at  rest  in  the 
gas  phase  has  zero  energy  by  definition.  Three  steps  s^,  and  are 
necessary  to  take  an  ion  from  the  gas  phase  to  the  surface  phase  as 


follows: 


S^:  As  shown  in  Figure  4,  an  energy  -ehd[N^]  is  required  to 
remove  the  ion  from  the  gas  to  a  point  remote  from  the  finite  dipole  layer 

of  adions.  In  Reference  2  a  discussion  is  presented  of  this  figure.. _ 

step  Sg!  The  energy  is  needed  next  to  move  the  ion  up  to  the 
bare  surface,  as  shown  in  Figure  2. 

Step  s=:  Finally,  there  is  the  interaction  energy  H  associated  with 
the  dipoles  as  the  particles  are  assembled  in  the  areas/.  This  energy  also 
includes  the  energy  of  depolarization.  ^ 


The  energy  of  Steps  S^  and  S*  for  Ions  is  N_^  <b 


si  f 

2  *v*  "si  "ei'io 


and  appears  in  the  exponentials  of  Equation  14.  The  energy  associated  With 

J  ,  N 
sa’  si 


S^  appears  in  the  configuration  integral  2[Nca,  NcS  3 


N  .  N 

si  sa 

(dax  )  ~r^r  ' j3 

n 


(d3x^)  exp(e 


...  (15) 


o=l 


n=l 


with 


...  (16) 


H  1  E  MoM  %  ,Ai) 

Hs  -  2  lSi<jSNsi  ion 

atoms 

where  r.  .  is  the  distance  between  particles  i  and  j,  M  and  M=M  -&£?  are  the 
ij  “  '  -  v  ?  p  -  o 

dipole  moments  at  zero  and  finite  coverage*  a  is  the  polarizability,  and  £ 
is  the  electric  field  given  by 

...  (17) 


C  =  E  M 
i^j  r*" 


in  the  dipole  approximation  to  be  used.  The  dipole  moments  have  been 
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replaced  by  their  average  values *  The  dipole-dipole  term  involves  MqM, 
rather  than  MM,  because  of  the  energy  involved  in  depolarization#  Both  the 
atoms  and  ions  have  been  assumed  to  have  identical  hard  core  potentials 
given  by  u(ij ) #  The  interaction  of  the  particles  with  the  metal  is  assumed 


to  be  uniform;  that  is,  no  lattice  sites  exist# 

The  configuration  integral  can  be  expressed  in  a  different  form**  which 
i «  amenabl e  to  approximate  calculations:  - — — 


€[N  ,N  .  ]  »</. 
sa  si  f 


N  .  +N  -N  .cp/2kT 
si  sa  si 

9 


...  (18) 


N_ 


/  s 

where  &  f  is  defined  by 

1  N 


.  N 


n^l 


s  -£u(ij) 

( da  x  )  e 


and  is  the  configuration  integral  for  the  hard  core  problem.  The  energy 
<p,  defined  by  Equation  18,  is  related  to  the  average  energy  E  , 

f  Nsa  Nsi  -H  /kT 

J  TT  7T  <**„>  V  5 

£  _  n=l  m=l  .  . . . . .  .  ,  ...  (20) 

6  f  Nsa  Nsi  -H  /kT 

4  7f  ( da  x  )  7f  (d  x  )e  B 
n=l  m=l 

by 

x  fN  .  <Pl 

*  It  [zkT-]  *  *9  ••• 

The  chemical  potentials  u  .  and  n  can  be  calculated  from  Equations  12,  13, 

SI  39 

14  and  18 


N 


si  A3 


Bs* 


if 


m  <b,  -  eAd  - 

lo 


n5iki  Ut 
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U  S  kT  In 
sa 


NsaAa 


-  (0  +el)  - 


a  f 


ae 


N  kT  M*  ,  > 

sa  f  .  1  o 

+i  s? 


sa 


sa 


The  departure 
22  and  23. 


rates  v.  and  v  are  then  derived  from  Equations  2,  3»  7,  8, 
1  & 


1  (N  ,cp) 

2  *N 


<*.  +eA#.i  "’si 


si 


„  .  “.iffl-  ,  I  !v»  x4  ' 

».  ■  r  1/a  ■ 77  “»  -37SC 


exp 


kT 


•“4  IT-  <»*♦>' 

_  _ sa 

kT 


In  Reference  1  the  f-factor  was  defined  by 

0.  s  0.  +  efdij  m 

i  io 

A  comparison  of  Equation  26a  with  the  last  factor  of  Equation  24  yields 


f  5 1  -  sfe*  (V*> 

The  temperature  dependence  of  0^  (B^  in  Langmuir  and  Taylor’s  notation)  is 

so  weak  for  Cs-[w]  that  it  was  not  observed  by  Langmuir  and  Tyalor.  The 

expression  of  Equations  24  and  25  in  the  form 

log  v.  =  A.  »  B,  /T 
-  i,a  i  ,  a  j.  |  a 

is  not  unique  when  B.  is  a  function  of  temperature.  In  this  case,  the 

i  >  a 

departure  rates  and  should  be  compared  with  experiment  rather  than 
the  non-unique  A  and  B  coefficients. 


A  different  expression  for  v  can  be  obtained  by  combining  Equations  10 

61 


and  11 


N 


jf-  a  t  exp  [e  (I- 
i 

Then  Equations  24  and  27  yield 
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with 


v 

a 


N 


2  f. 


si 


exp 


N 


si 


3r  55- 


a*  \ 

exp 

r»T 

i  ' 

i  / 

0  =  0.  -  el  +  60 

a  i 

4 

Langmuir  and  Taylor  have  used  forms  such  as  Equations  24  and  28  as  indi¬ 
cated  fey  the  fact  that  their  B_  and  B.,  differ  by  ln2, 

'  OL  as 

The  ratio  of  the  number  of  adatoms  to  adions  can  fee  obtained  by  com¬ 


paring  Equations  25  and  28 


with 


w 


6. 

i 


N  ,/N  , 

si  so’ 


N  /N 
sa  so 


...  (29) 


N  being  the  number  of  cesium  particles  in  a  monolayer.  This  equation  can 
60 

also  be  obtained  directly  from  the  condition  for  chemical  equilibrium, 


“sa  =  *si  +  V 

BRAGG -WILLI AMS  APPROXIMATION 

An  approximate  form  for  ,  in  Equation  18  can  be  obtained  with 

sa  si  . 

a  Bragg-Williams-type  approximation. *  First,  the  factor  w  ^  is  evaluated 
for  low  surface  densities: 


0  N  +N  .  i»  « 

<7  >  **  st  =  J  w 


N  .  N 

81  (d3x  )  tt 

-  ffi  - 

m=l  ~  n=l 


3~  '  -sa  (d3x  )  e 


-Eu(ij)/kT 


N  *N  . 
sa  si 


1  Nrrr*2,  Nsa+Nsi 


0 


...  (50) 


with  r  being  the  hard-core  radius  of  u(ij).  Thus 

1 

2 


=c/-iNnr*2 
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The  r*  was  chosen  as  in  Reference  2  So  that  <?i nr*3  a  2,  where  =  density 
of  a  monolayer  for  Cs-[w]  a  3.563  x  lO1**  cm“2,  and  then  r*  a  4.2  x  10“®  cm. 
Consequently , 

<2/  =  ^(i-0)=  «/(i-e.-e  )  ...  (31) 

1  -  1  a 

and  the  "entropy’'  factor  in  the  departure  rates  will  have  the  same  form 
that  Langmuir  and  Tyalor  obtained,  that  is,  the  A.  and  A  factors.  Equation 

1  Si 

31  may  he  compared  with  the  more  accurate  expression  of  Tonks^ 

/  1  1.30763  +  0. 30?Q4 

**  t  ~  1  ♦  1.8146  +.  2.5736s 


Equation  31  is  assumed  to  be  true  for  all  coverages,  although  this  is 
rigorously  true  only  for  small  coverages. 


The  energy  cp  calculated  in  Equation  18  is  calculated  from  the  approxi¬ 
mate  form  of  E  obtained  by  the  firagg-Williams  approximation .  E  is  first 

8  S 


rewritten 


E  =  1  N . 

S  2  S3 

0 

L 

...  (32) 

tt 

11 

i& 

g.(Xx  ,  Xa ) 

-  ..  x  »8i 

,  xa  )  -j— 

...  (33) 

g  ,Cx  i  ,  Xa  ) 

1  _ 

a_  -HrM 

xNe  s 

...  (34) 

where  EL(xi ,  xa )  is  the  interaction  energy  of  ions  1  and  2  according  to 
Equation  16.  The  Bragg-Williams  approximation  consists  of  setting  the  pair 
correlation  function  g;(xi  ,  xg )  in  Equation  33  equal  to  that  for  a  random 
distribution  of  particles  outside  the  hard-core  radius  r*.  One  then  ob¬ 
tains  from  Equations  21,  32,  33  and  34 


N 


si  1 


M  M 
o 


2  P" 


2  nrdr 


...  (35) 


AI-7979 

3.7 


The  electric  field  &  for  this  distribution  is  given  by 


and  then 


M  & 


M  a  M  - 
o 


N  2ttN 

-4*  2nrdr  a  ,-f1  M 


71 


M 


2ttN 


l  +  o 


si 


The  chemical  potentials  can-  now  be  calculated  from  Equations  22  and  2$' 
with  Equations  31 »  35  and  37 


N  .A3  / 

M> a  kT  in  -  d.  -  eA*  +  J™  0.kT  * 

si  q^ef  ^  io  •'f  i 


tp 


**fP«*'i 


N  As  / 

U  a  kT  In  r^-j-  -  («J  +  el)  +  £-  6  kT 

sa  da  of  *  ao  <*  *  a 

a  i  i 


Similarly,  the  departure  rates  are  Obtained  from  Equations  24  and  25 

N 


vi  •  ft  exp(^  ei  exp  /  *io *  eA®  ‘ 


JL 


kT 


7 


_  sa 

Va  *  fa  W~  exP 


*ao 

**?  I  *  "kT 


...  (41) 


This  work  function  change  A -5  is  given  by 


2na  M  9.  4  M  9 .  __ 

A 0  =  *  - 2  - ;  «  sa  2221 

P*a  ■>  ■  r 


1+  <*9.  r*s  1+  o-e. 


(42) 


It  will  be  shown  that  9^0  for  0<O.5O.  The  constants  O’  and  4M  /r*3  in 

—  o' 

Equation  42  have  been  determined  previously  in  Reference  2  to  be 
4M 


pr  =  10.8  ;  o  =  2 


AI-7979 

l8 


The  f- factor  is  obtained  from  Equations  26b  and  35 


f  m  1  - 


SN  . 
si 


=1-0.4 


1+-  9 
<P  1  2  i 


1+50. 

1+0. 
_ i 

rJis" 


=  i  - 


M 

o 

2er* 


1+0, 


**.  (43) 


Thus  the  f-factor  increases  slowly  from  0.60  at  =  0  to  0.73  at  ©^  =  1*00. 

n  _  . 

Experimentally, '  f  decreases  from  0,8  at  ©  =  0  to  0.65  at  6  =  0.80,  For 


Cs-Cw],  6-0^  for  6<0.55»  The  above  discrepancy  is  considered  minor  since 
the  main  dependence  of  ^  on  9  occurs  through  the  factor  44  in  Equation  26a, 
f  being  approximately  constant. 


CALCULATION  OF  40(9 ,T) 


This  section  is  divided  into  two  parts.  First,  the  modifications  of 
1  2 

the  earlier  theory  *  in  the  light  of  the  theory  just  developed  are 
described.  In  the  second  part  the  effects  of  adatom  polarizability  are 


investigated. 

Modifications  of  Earlier  Theory 

The  work  function  Change  is  given  implicity  by  the  two  Equations  29 
and  42  in  the  two  unknowns  ©^  and  40  for  given  9  and  T.  Previously,  in 
Reference  1,  the  atom  to  ion  ratio  was  assumed  to  be  given  by  a  simple 


Boltzmann  factor  with  a  degeneracy  factor  of  2.  The  correctness  of  the 
factor  2  is  very  questionable,  according  to  Equation  29,  based  on  a  chem¬ 
ical  equilibrium  for  Cs**Cs+  +  e“,  the  degeneracy  factor  is  given  by 

</e. 

1 
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Before  discussing  Equation  29  further,  it  is  interesting  to  evaluate 
as  a  function  of  0  from  Langmuir  and  Taylor's  experimental  A0(0)  with 
Equation  42 

10.8  0. 

40  =  volts  or  0±  =  -fgy-. 


The  results  are  shown  in  Table  1. 

TABLE  I.  ADION  COVERAGE  9  VS.  TOTAL  COVERAGE  0 _ 

FOR  ZERO  ADATOM  POLARIZABILITY 


A  0(6) 

=  0.8?4 

IH! 

. 

vjv 

O 

O 

2.003 

2.420 

2.73 

2.83  ; 

1*89 

2.92 

2.92  ; 

’*90 

2.87 

0 

s  .10 

*20 

.30 

.40 

.50 

.55 

.60 

.65 

*70 

.75 

.80 

e. 

i 

=  .097 

.192 

.295 

.406 

.511 

.551 

.576 

.589 

.589 

.580 

.567 

0/0 

i' 

0,-0 

=  .975 

.96 

.98 

1.02 

1.02 

1.00 

*96 

*91 

.87 

*77 

.71 

i  a 

=  .11 

.25 

.43 

.67 

1.00 

1.22  3 

L.38 

1*51 

1.59  3 

L.64 

1.67 

The  difference 

between  0. 

and  0 

for  0< 

0.55  is 

not 

considered  to 

be  caused 

by  adatoms;  for  example,  patchiness  of  the  tungsten  filament  might  be  im- 

n 

portant,  but  has  not  been  considered.  The  strong  decrease  of  0^/0  for 
0>O.55  is  considered  to  be  real.  It  might  be  noted  now,  however,  that  in 
the  second  part  of  this  section  the  inclusion  of  adatom  polarizability 
leads  to  many  less  adatoms.  In  both  cases,  the  rapid  turnover  of  A0(6)  for 
0X),5$  is  to  be  attributed  to  adatoms. 

In  Equation  29,  at  the  lower  coverages  the  energy  factor  I  is  domi¬ 
nant,  and  6  <<6  ,  At  the  higher  coverages,  0>O.55,  the  preceding  table 
shows  that  (0.-0  )/(l-0)  is  slowly  varying,  but  so  is  E  =  E  -(l-f )Ad  since 
Ad  is  near  maximum.  At  these  higher  coverages,  the  model  is  uncertain 
because  of  other  factors  such  as  lateral  bonding  of  cesium  so  that  the  des¬ 
cription  of  simple  dipoles  may  no  longer  be  valid.  It  is  therefore 
practical  to  take 
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N 


-  A  exp  (- 
si 


E 

kT 


(44) 


With  A  considered  a  constant  independent  of  substrate  material. 


In  the  previous  analysis,  the  electric  field  was  taken  to  be 


£  = 


(l-f)40 

er. 

i 


..  (45) 


which  is  to  be  compared  with  Equation  36.  Equation  11  Of  Reference  1 
results  if  the  electric  field  is  given  by  Equation  45  rather  than  Equation 
36.  In  solving  Equation  21  of  Reference  1  for  40(9 ,T ) ,  the  additional 
6-dependence  through  the  (l-f )  factor  is  not  necessary  but  does  help  in 
obtaining  a  maximum  in  40(9) .  The  (l-f)  factor  (taken  to  be  that  experi¬ 
mentally  determined)  in  the  energy  E  of  Equation  44  is  actually  sufficient 


to  obtain  the  40 -maximum.  The  (l-f)  factor  in  the  electric  field  in 
Equation  45  is  not  on  firm  theoretical  ground  since  the  term  f40  can  be 
considered  to  consist  of  two  parts:  1)  the  energy  to  rearrange  the  surface 
particles  from  the  initial  to  the  final  configuration;  and  2)  the  energy  to 
bring  the  ion  to  the  surface  under  the  electric  field  of  the  final  configu¬ 
ration.  Despite  the  above  uncertainties,  the  electric  field  expression, 
Equation  45,  will  be  maintained  since  the  (l-f)  factor  is  varying  slowly 
and  any  resulting  changes  in  40(9, T)  from  the  use  of  Equation  45  instead  Of 
Equation  j6  should  be  minor.  Two  effects,  atomic  polarizability  (to  be 


discussed  below)  and  lateral  bonding,  have  been  neglected  so  that  the 


40-results  at  high  coverages  will  be  only  semi-quantitative .  The  results 


presented  in  Reference  1  must  now  be  reinterpreted  slightly  with  the 
following  ambiguity.  From  Equation  44  and  the  expression  obtained  in 


Reference  1  for  N  /N  .  which  fits  the  Cs-[w] 
"  -  sa'  si 


N 

N 


sa 


=  Aexpl- 


Vu" 


Si 


kT 


-  2  exp 


l.O54-(l-f)40\ 

,052 


(46) 
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and  thus  with  A  -  2e“  it  is  necessary  that 
*  +  19.2  Eq  =  20.3 

for  example  with  =  0,  E„  =  1.054  ev.  For  =  1,  Eq  =  1.01  ev  and  for 

'i'  =  -1 ,  E  a  1.11  ev.  Thus  for  A  not  differing  from  2  by  orders  of  magni- 
o 

tude,  Eq  remains  approximately  1  ev.  In  this  manner,  the  rapid  turnover  of 

A0(0,T)  with  respect  to  6  at  approximately  50%  coverage  can  be  explained  by 

the  existence  of  adatoms  with  an  atomic  energy  of  evaporation  - 

0  ■  0  -1  +  0.  -E  a  4.62  -  3.87  +  2.04  -  E 

ao  o  xo  o  o 


a  2.79  -  E 

o 

8 

which  is  in  fair  agreement  with  deBoer's  rough  estimate,  0  1.4  ev. 

£LO 

Inclusion  of  Atomic  Polarizabilities 


It  is  to  be  expected  that  the  adatoms  will  be  polarized  by  the  electric 
field;  therefore,  the  adatoms  will  decrease  directly  the  work  function 
change  40.  This  effect  is  easily  included  in  the  theory ,  and  the  basic 
equations  now  become 

...  (47) 
...  (48) 
...  (49) 


M.  =  M  -  a£ 
X  o 

M  =  •  a  € 
a  a 


Ad  a  2WJJ  (M.0.  +  2M  0  ) 
*  x  i  a  a 


2m, 


<MA  *  w 


...  (5P’ 


^  .  Ae-E/KT 

i 


in  the  Bragg-Williams-type  approximation ,  where  o  is  the  atomic  polariz- 

a 

ability  and  the  factor  2  arises  in  Equation  49  because  the  dipole  lies  out¬ 


side  the  metal  surface.  At  zero  coverage ,  the  adatom  is  assumed  to  have  a 

negligible  dipole  moment.  Without  the  Bragg-Williams-type  approximation , 

the  induced  dipole  moment  M  will  depend  on  the  location  of  the  atom 

a 
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At  high  coverages  (0>O.55)  where  we  shall  find 


relative  to  the  ions. 


adatoms  present,  each  adatom  has  two  or  more  adions  as  nearest  neighbors, 
The  surface  is  sufficiently  crowded  that  the  electric  field  in  Equation  $0 
is  a  good  approximation.  At  low  coverages  adatoms  might  fee  expected  to  fee 
found  adjacent  to  adions  according  to  the  attractive  dipole-dipole  inter¬ 
action  (in  Reference  8  a  different  reason  is  given) .  In  this  event,  it 
should  be  noted  that  by  electrostatics  the  electric  field  at  the  center  of 
the  atom  due  to  the  adion  and  its  image  has  its  vertical  component  pointing 
towards  the  substrate  rather  than  away  as  depicted  by  deBoer.  deBoer  has 
recently  discarded  this  picture  of  adatoms  adjacent  to  adions  and  now  con¬ 
siders  the  adsorbed  cesium  as  one  species  whose  adsorption  bond  changes 

9 

with  coverage. 


Equations  47-50  can  be  combined  to  obtain 

A0  =  2to1M  0  1~°raC9 

°  1  l+QrC0.+a  C& 
l  a  a 


which  should  be  compared  with  Equation  42.  Equation  52  offers  another 

means  for  obtaining  the  AO  maximum.  Equation  52  can  be  used  to  obtain  0^, 

and  0  =  0-0. ,  from  the  experimental  A0(0)  if  the  adatom  and  adion  polariz* 
a  x 


abilities  are  assumed  equal.  Table  II  was  obtained  with 

i-e+0. 

"  •  lo-80i  ml 


(52) 
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TABLE  II*  ADION  COVERAGE  6.  VS.  TOTAL  COVERAGE  6  FOR 

■  -  -  -  i  • 

EQUAL  ADATOM  AND  ADI ON  POLARIZABILITIES 

9  0.10  0.20  0.30  0.40  0.50  0.55  0.60  0.65  0.70  0.75  0.80 

Ad  .8738  1.500  2.003  2.420  2.73  2.83  2.89  2.92  2.92  2.90  2.87 

0±  .0974  ,194  ,297  *402  *504  *550  .593  *833  *669  *704  *737 

ei/e  *974  ,970  .990  I. 01  1.01  1.00  .988  ,974  ,956  ,939  *921 

The  Ad-maximum  is  obtained  with  less  adatoms  than  in  fable  I*  the  reason 

being  that  the  adatoms  are  directly  decreasing  Ad  with  their  induced  dipole 

moments.  Previously,  the  adatoms  caused  the  decrease  of  Ad  because  of  the 

(l-f )  factor  in  1  in  Equation  46.  The  values  of  9^  are  consistent  with 


i 


'  E 

=  2  exp 

’  kf 

=  2  exp 

l  i 

1.19  ev  -  (l-f)Ad 
*052  ev 


where  (l-f )Ad  is  taken  from  the  experimental  data*  This  form  for  E  is  not 
really  correct  since  no  energies  involving  the  adatom  dipoles  are  included* 
These  interactions  can  be  included  by  rederiving  Equations  33  and  35  to 


Obtain  a  new  expression  for  m 


d8r  gftCr)  Hft(r) 


(54) 


where  gi(r)  is  the  ion-ion  pair  correlation  function  given  by  Equation  34, 
and  g  (r)  and  g  . (r)  are  the  atom-atom  and  atom-ion  pair  correlation 

®  AX 

4 

functions  given  by  expressions  similar  to  Equation  34.  EL  (r)  is  the  same 

as  in  Equation  34.  H  (r)  and  H  . (r)  are  the  atom-atom  and  the  atom-ion 

a  . ai  ~  ‘  '  — 

interaction  energies,  respectively*  Equations  18,  21-26  and  32,  which 
involve  and  5>,  still  hold.  The  Bragg-Wiiiiams  approximation  for  the  pair 
correlation  functions  could  be  made  in  Equation  54  and  the  evaporation 
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rates  v  ,  v.  and  v  calculated.  These  calculations  have  not  yet  been  made, 
a  i  e 

but  would  be  relatively  simple.  It  is  very  likely  that  Other  forces 
besides  the  dipole-dipole  forces  envisioned  in  Equation  3  b  will  be  very 
important  at  high  coverages  such  as  cohesive  metallic  bonds,  the  calcu¬ 
lation  preceding  from  Equation  3b  would  then  net  have  much  significance . 

The  significant  factor  shown  in  this  section  is  the  possible  effect  of 
adatom  polarisation  on  the  work  function.  It  must  be  remembered  that  the — 
existence  of  adatoms  and  adions  as  two  distinct  species  is  a  debatable 
point.  It  might  be  expected ,  however,  that  the  vibrational  spectrum  of 
adatoms  and  adions  would  be  coverage®  and  temperature-dependent  and  that 
these  effects  could  be  observed  in  infra-red  adsorption  or  reflection 
Spectra.  Adsorption  as  two  species,  then,  would  not  have  to  remain  a 
speculative  point. 

SUMMARY 

12  3 

The  previous  ONE  work  *  might  be  summarized  as  follows.  In  Refer¬ 
ence  2  a  Statistical  mechanical  treatment  yielded  evaporation  rates 
and  in  terms  of  various  microscopic  constants.  Examination  of  these 
constants  showed  how  to  treat  different  substrate  materials.  The  treat- 
ment  yielded  the  same  form  for  the  A-factors  (v  =  eAe*  '  )  that  Langmuir 
predicted  (the  basic  assumptions  were  the  same,  but  the  treatments  were 
quite  different  in  the  two  oases).  The  energy  factor  B  was  obtained  from 
outside  the  statistical  mechanical  formalism,  whereas  Langmuir  did  not 
derive  an  expression  for  B.  No  adatoms  were  assumed.  In  Reference  1 
adatoms  were  found  capable  of  explaining  the  maximum  in  the  work  function 
change.  By  this  means,  predictions  of  d<S(0,T)  were  for  arbitrary  substrate. 
Theoretical  expressions  for  va  and  were  calculated  from  the  Cs-fw]  data 
by  assuming  that  a  variation  of  the  substrate  work  function  left  Aft,  A^, 
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unchanged  and  changed  in  a  manner  indicated  in  Reference  2.  In  Ref¬ 
erence  3  a  calculation  of  the  f-factor  0^  =  0^  +  fA0)  was  made*  This 

calculation  did  not  include  any  energy  associated  with  the  rearrangement  of 
the  addons  as  the  test  ion  is  brought  to  the  surface*  The  results  did 
agree  well  with  the  experiment* 

The  present  treatment  calculates  the  expressions  for  the  evaporation 

rates  v  ,  v.  and  v  with  adatom  and  adion  adsorption.  For  polarizable  _ 

adatoms*  the  calculation  of  v  and  v.  is  only  indicated*  In  this  treat- 
ment ,  a  general  expression  for  the  f  factor  falls  directly  out  of  the 
formalism,  in  the  Bragg-Williams  type  approximation,  f  (6^)  does  net  agree 
with  experiment'  as  well  as  f (9 )  in  Reference  3  did,  but  the  agreement  was 
still  considered  satisfactory*  The  expressions  for  and  differ  from 
those  obtained  by  Langmuir  for  appreciable  adatom  coverage.  This  differ¬ 
ence  is  not  in  the  right  direction  to  explain  the  small  4.8(9-0s/2)  depend- 

4 

ence  found  experimentally  in  A  in  addition  to  that  prediction.  Finally, 

El 

it  was  found  possible  to  explain  the  maximum  in  A0  with  adatom  polariza¬ 
bilities  whereas  previously  in  Reference  1  it  was  necessary  that  the 
f  factor  decrease  with  8. 

There  are  many  points  at  which  the  present  theory  may  fail  and  these 
have  often  been  pointed  out  in  the  text.  One  point  often  made  was  that 
lateral  bonding  of  the  cesium  and  possible  accompanying  effects,  such  as 
change  of  polarizability,  could  greatly  affect  the  calculations.  Patchi¬ 


ness  of  the  surface  was  never  included. 

Most  useful  would  be  measurements  of  v  ,  v  and  v 

-  -  -  ~  -  -  a  i  e 

and  T  for  a  low  work  function  substrate  material. 


as  functions  of  9 
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IV.  SPACE  CHARGE  ANALYSIS  FOR  LOW  PRESSURE  THERMIONIC  DIODES 

C.  Warner,  III 

INTRODUCTION 

Space  charge  analyses  of  collisionless ,  plane  electrode ,  thermionic 

1=8 

converters  have  been  performed  by  various  authors.  These  analyses  can 

9  . 

be  considered  as  extensions  of  the  work  of  Langmuir  who  solved  the  space 

charge  problem  with  one  charged  species.  If  both  ions  and  electrons  are 

emitted ,  the  potential  distribution  can  assume  a  variety  of  forms.  This 

paper  first  considers  the  space  charge  problem  for  infinite  spacings. 

Particular  types  of  solutions  are  shown  to  exist  in  different  regions  of 

fl  .  P  .  space  as  shown  in  Figure  1,  in  which  P_  is  the  dimensionless 
m  min  m 

potential  where  the  electric  field  first  vanishes  (this  point  is  sometimes 

an  inflection  point) :  and  Tl  .  is  the  potential  of  the  collector.  These 

mm 

potentials  P  and  P  .  .  are  shown  in  Figure  2.  In  Figure  1  it  Is  shown  that 
‘  m  min 

the  independent  paramaters  and  "^min  determine  the  form  Of  the  potential 

distribution.  The  parameters'!)  and  P  ,  also  determine  the  volt-ampere 

m  min 

curves  for  given  emission  conditions  (that  is*  where  8  is  constant) .  Each 
point  on  the  volt-ampere  curve  can  be  identified  with  a  particular  form  of 
potential  distribution,  and  thus  the  continuous  change  of  potential  form 
becomes  apparent  as  the  collector  potential  P^  varies.  Finally,  the 
effect  of  finite  spacing  is  discussed. 

THE  BASIC  EQUATIONS 

The  potential  distributions  are  obtained  by  solving  Poisson's  equation 
in  which  the  charge  density  is  given  in  terms  of  the  potential.  This  ex¬ 
pression  for  the  charge  density  depends  on  the  particular  potential  shape 
chosen.  The  motion  of  the  charged  particles  is  assumed  to  be  governed 
entirely  by  the  potential  distribution  P,  that  is,  the  Vlasov  equation  is 
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According  to  the  Shape  of  the  Potential  Distribution 
(The  Numbers  1  to  7  Correspond  to  Curves  1  to  7  in1  Figure 
“H  and  1  •  are  defined  in  Figure  2) 


CURVE 


Al-7979 

30 


Figure  2.  Some  Typical  Potential  Distributions 
for  Infinite  Spacing 


being  solved*  The  emitted  charged  particles  are  assumed  to  have  a  Max¬ 
wellian  distribution  corresponding  to  the  emitter  temperature* 

For  the  potential  distribution  shown  in  Figure  3*  Poisson's  equation 
taxes  the  form 

^  _  _ _  _ 

a  {  1-erf  v/iT  +  erf  }  -  ^5-  |1  ±  erf  ^Tln-1)  J 

..*  (1) 

the  plus  sign  holding  for  o<§<?  and  the  minus  sign  for  S  <5<§  .  The  dim- 

m  no 

ensionless  potential  is  defined  by  1)  s  eV/kT ,  and  vanishes  at  the  emitter, 

V  being  the  potential  and  T  the  temperature*  The  parameter  8  =  n  Q/n+0 

where  a  /2  and  n  /2  are  the  emitted  electron  and  ion  densities  at  the 
-o  +0 

emitter.  The  dimensionless  distance  is  measured  in  ion  Debye  lengths 
4tfn  e3 

I?  _  x3  x  =  distance  ...  (2) 

Equation  1  is  typical  of  the  form  of  Poisson's  equation.  In  all  Cases,  the 
first  integration  may  be  done  analytically.  For  the  potential  distribution 
of  Figure  3*  there  results 


where 


|ff,  F,  (T)  -  f,  <y 

*  h  W  *  Fa  <y 

«  F,  (T)  -  F_  U\) 
P  P  m 

*  F4  (T!)  -  f4  <-y 


0<?<§ 

...  (3) 

m 

s  <s<s 

...  (4) 

m  0 

!<!<!„ 

P  ffl 

...  (5) 

t  <§<  5 

m  max 

...  (6) 

Fx  (T) 

f2  (d) 

F3  (H) 


-T|  -T) 

2e~  -e  “  G(T)  -H)  ♦  200(11)  -  8e  0  Q(T|  +H) 

m  m 

...  (7) 
-n 

e  m  G(T1  -11)  +  206(11)  -  8e  “  6(1)  +D)  ...  (8) 

m  ~  '  _  ■ 

-n  - 

e  m  6(11-11)  *  20  e'  -  0e  “  6(Hm+H)  ...  (9) 

n  n 


Ai-7979 

31 


1 


I 


-T]  -V.  _ 

F.  (1i)  =  e  m  G(\*U)  +  Be  m  G(\+T|)  ...  (10) 

«  m  in 

G(T|)  S  e^  erfc  +  £_  */n""  ...  (11) 

v'rt 


Exhibiting  Potential  Maxima  »  ^m  and  , 

and  Potential  Minimum  T|  .  max 

min 

The  first  integrations  for  all  the  various  potential  forms  are  basically 
contained  in  Equations  7  to  10.  For  example,  the  potential  distribution 
for  curve  1  in  Figure  2  is  given  by 

( 3r)  *  (d?)  „  ,  *  *1  (1'\  *  W\  *  0) 

■  min 

*1  ^~\in  ^“^min  T]  =o)  ...  (12) 

ffl 

which  is  obtained  from  Equations  3  and  7.  As  discussed  later,  the  deriv¬ 
ative  (dT|/d|)9  ■p_>n  must  be  chosen  properly  to  obtain  the  vanishing  of 

|  min 

(dT!/d? )s  [1,T1  ,0]  at  T1=T|  . 

min  m 

The  parameters  B ,  T|  ,  T)  and  H  are  not  independent  for  the  distri- 
*  m  m  max  -  -  . . 

bution  of  Figure  3?  The  continuity  of  the  electric  field  implies 

Fj(0)  .  F2(\)  •  Fj(0)  .  F?(.\)  ...  (13) 
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from  which  is  obtained  8=9  [11  ,T!  ].  The  vanishing  of  the  electric  field  at 

jn  m 


max 


F,(n  )  -  F,  (-T1  )  =  o 
4  max  4  m 


from  which  is  obtained  T1  ...  =  T|  !!S,H  71  ],  Similar  relationships  arise 

max  max  ’  m’  m  ~  * 


is  obtained  by  reflection  of  the  Solution  between  5 
max  m 


for  Other  potential  distributions.  It  should  also  be  noted  that  the  so* 
lution  far 

max 

and  t  .  about  §=?._  ,  2§ .  -  I  ,  etc. 

_  max  max-  max  - 

SOLUTIONS  FOR  INFINITE  SPACING 


The  oscillatory  solution  shown  in  Figure  3  corresponds  to  infinite 
spacing.  In  Figure  2  are  shown  other  possible  forms  for  infinite  spacing*. 
Solutions  2-5  are  also  solutions  for  finite  spacing  if  the  collector  is 
placed  at  any  point  at  which  all  the  ions  and  electrons  are  moving  towards 
the  collector.  Besides  the  oscillatory  solutions,  potential  distributions 
such  as  curves  1,  6  and  7  in  Figure  2  are  possible*  The  break  in  curves  1, 
6  and  7  corresponds  to  an  infinite  distance*  To  obtain  these  solutions, 


the  input  parameters  must  be  selected  so  that  the  second  integral  to  obtain 

(dT|  V® 

dflH  fl 


in 


'min 


Equation  12  is  taken  equal  to  the  negative  of  the  minimum  of  F^CT],^  = 

T!  .  ,  T|  =  0)  -  F.  (T'sTI  .  ,T)  =  T|  f  Tl  =0)  with  respect  to  H.  The  elec« 
min’  m  1  min’  m  min’  m  - 

trie  field  d’W/dl  vanishes  at  the  corresponding  "plasma  potential"  and  the 

second  integral  for  I  diverges  logarithmically*  A  proper  choice  of  8  for 

given  values  of  T1  and  T)  yields  curves  6  and  7.  That  is,  with  proper 
m  m 

/dT]  3 

choice  of  8,  exhibits  a  minimum  with  respect  to  11  and  also  vanishes  at 

that  minimum  point*  The  second  integral  for  f  also  diverges* 

Systemization  of  Space  Charge  Solutions 

The  various  types  of  solutions  shown  in  Figure  2  are  used  to  define 

regions  on  the  H  ,  T|  .  plot  in  Figure  1.  Lines  of  constant  8  are  also 
m  min  ~ 
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shown.  The  lines  labeled  with  numbers  2 ,  5  and  7  and  the  regions  with  the 

circled  numbers  1,  3»  4  and  6  can  be  identified  with  the  curves  in  Figure  2 

with  the  same  numbers.  The  B— values  in  Figure  1  are  given  by  the  following 

For  0<Tlmin  <'Flm  (e.g.  curve  3  in  Figure  2) 

G(T1  .1  .  )-i 
a  m _ m  mm 

=  GCH  )-G(T  ) 

m  mn 

For  ('ll  )  <T  =  -  T!  so  (e.g,  curves  4  and  5  in  Figure  2 

max  in  min  m  ° 


Be 


QCH  +n  )-i 
m  m _ 


...  (16) 


2G(T  )-e  m  [*G(T1  +\)+l 

m  mm  - 


The  function  Tlmny(Tlm)  is  defined  below.  Equation  16  is  equivalent  to 


Equation  13 


For  -  oo<T|  *  „  =  -T:  <11 _ (e.g.  curve  6  in  Figure  2) 

min  m  max 


Be 


1L 


m 


G(tl  -T)_i 
in 


s  min 


2fG(t  )*e^l  +  e  “  rG(T|  +T1)-G(T1  +T,  )] 

m  m  m  m 

The  minimum  is  with  respect  to  t|  for  ■»'fl  <tKO*  Equation  17  is  derived  from 

m 

-H  r  -T) 


...  (17) 


($[  *  •"'*  - » [♦’ "  {«vil)**<Vv} 

+  2CG(H  )  -e  ]  TKO 
m  -  4 

*  M1L,H)  .  .  b|4  .  « 


...  (18) 


m*  m ’  m ' 

The  function  A/B  has  a  minimum  with  respect  to  ”11  and  Equation  17  therefore 

leads  to  the  divergence  in  the  second  integration  for  infinite  spacing. 

For  \(\j_H\<\,.,_(e.g.  curve  1  in  Figure  2) 
m  min  m  min  -  - 


Be 


2  -  e  ■»  m  g'(T  4  -T)  ) 

.  _  _  , ,  ,  ~  ,  min  _k 

G-cn  ) 

m 
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In  this  region,  (dT]/d§ )"  is  given  by  Equation  12.  According  to  the  defi¬ 
nition  of  as  given  in  the  introduction  and  in  Figure  2,  Equation  19  is 
just  equivalent  to 


SF 


(*m=t|  ,  1  =  H  .  ,  T)  =  0)  =  0 

m’  m  min’  m 


It  is  also  of  interest  to  give  the  equations  for  the  boundary  lines. 


For  T|  =  T)  (T>  )  (i.e.  curve  8  in  Figure  1  or  21. 

m  m  min 


The  function  )  is  given  implicitly  by 


He  m-l)  -  e  miB  m  G(H  .  )- 

_ _ _ _  _ tun 


lmin“  ja? 


G<*n  )_i 
m 

#(1. ,  -ft*)  . 

°'<V 


which  is  obtained  by  eliminating  B  from  Equation  19  and  the  equation  that 

(dtv)3 


states  that 


=  0  at  the  emitter: 


F  (H  ,11  . .  ,8)  =  F,(0,T  .  ,S) 
1  m’  min*  1  *  min’ 


For  T|  =  Tl7  ( T:  )  for  H  =0  (See  Figure  3) 
m  mm  max 


The  curve  1)  =  Tl  (11  )  for  Tl  s  0  is  given  by 

m  mm  max 


m  m 


-r 


,)-e 


m 


G(T!  +11  )+l 
m  m 


m  m  m  m 

G(H')-G(0) 

m 


...  (21; 


which  is  obtained  by  setting  11  s  0  in  Equation  14  and  eliminating  B 

max 

between  Equations  13  and  14. 

For  =  H'  (H  )  for  T|.  .  s  11  (see  curve  5  in  Figures  1  and  2) 

min  max  m  min  max 


X  )-l 

m  max.  . 

-  -V  ) 

m  max 


»tl 


=  e 


max 


)  _  G(T|  _T  )-l 
in  m  max 


which  is  obtained  by  letting  approach  -^max  ih  Equations  15  and  14  and 
eliminating  S , 


The  solutions  for  TL._  =  0  are  not  shown  in  Figure  2,  but  serve  as  a 
good  starting  point  for  explaining  Figure  1.  The  discussion  breaks  natu¬ 
rally  into  the  three  regions ,  0<8<0.405,  0.405  <B<0. ?86  and  0. 786  <B<1 ,000. 

The  3=0.405  line  is  asymptotic  to  1i_=0  line  for  negative  The 

m  min 

B  =  (0.786)“^  line  is  tangent  to  the  11  sst]  .  line  at  D .  .  —.20  and  does  not 

m  min  min 

cross  that  line. 


The  Region  0.405<B<0.786 

The  B  a  0.5  line  is  considered.  As  the  collector  Voltage  11...  _  is 

nn 

varied,  the  potential  distributions  will  progress  through  the  sequence 

shown  in  Figure  2.  As  11  .  is  increased  from  zero  the  oscillating  region 

min 

in  Curve  5  of  Figure  2  decreases  in  amplitude  and  wavelength,  until  finally 

curve  2  is  reached  where  T  =T  .  .  As  1  .  is  further  increased,  the  poten- 

tial  distribution  resembles  curve  1.  For  B  =  0.5,T  varies  little  as  T  . 

m  “  Din 

goes  to  infinity.  If  ^  is  made  negative,  the  distribution  resembles 

curve  4,  of  Figure  2.  As  11  .  is  made  more  negative,  T  decreases  to 

zero  (T  =0  in  Figure  3)  and  finally  11  =11  .  (curve  5  in  Figures  1  and 

max  6  max  min  0 

2) .  Curve  5  has  E=d11/dS=0  at  a  finite  distance  and  everywhere  thereafter. 
Further  decrease  of  H  .  develops  a  collector  sheath.  As  T|  .  is  decreased 
T  finally  vanishes  and  curve  7  in  Figure  2  is  obtained.  To  proceed 
further,  the  following  symmetry  principle  must  be  recognized? 

Every  ion-rich  solution,  11=f(§),  with  8=8^<lt  cor¬ 
responds  to  an  electron-rich  solution,  H=»f(l),  with 
B=i/°^.  The  electron-rich  solution  is  the  negative 
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of  the  ion-rich  solution.  The  spacing  ?  in  the 
electron-rich  solution  is  measured  in  electron 


Debye  lengths  instead  of  ion  Debye  lengths  for 
the  ion-rich  solution.  Only  the  electric  charges 
of  the  ions  and  electrons ,  and  not  the  masses , 
enter  the  steady-state  equations. 


Further  decrease  of  for  8=0 . 5  thus  corresponds  to  an  increase  of 

ILi-  from  the  E(emitter)  =  0  curve  (curve  8)  in  the  first  quadrant  along 
min 

the  8=2  curve.  In  this  manner,  as  the  collector  sheath  voltage  increases 
without  limit,  the  "plasma  potential"  («-TI_  in  first  quadrant)  approaches  a 
finite  value  asymptotically.  In  accordance  with  the  symmetry  principle, 


Equation  17  evaluated  for  T)  =0  becomes  Equation  19  for  I)  under  the 

n  ni  in 

transformation:  TK-TL ,  8-*l/8  and  THN  .. 

m  min 


The  Region  0<4<0.405 

For  8<0.4O5,  as  ^min  goes  to  negative  infinity,  a  finite  barrier 
remains.  This  possibility  is  unexpected  since  even  with  electrons  only, 
the  barrier  can  be  made  to  vanish  for  any  finite  spacing.  This  infinite 
spacing  and  infinite  voltage  situation  is  perhaps  unphysical,  but  is 
resolved  by  reference  to  Figure  4a  in  the  section  "Emission  Requirements 
for  Removal  of  Space  Charge  Barriers,"  If  the  limit  a>  is  taken  first, 
then  the  limit  results  in  zero  space  charge  barrier  only  if  8<0,4049. 

If  oo  first,  then  the  space  charge  barrier  is  always  removed.  It  must 
therefore  be  remember  that  the  formation  of  a  "plasma"  is  not  the  only 
means  by  which  an  infinite  spacing  solution  may  be  obtained.  The  Other 
means  corresponds  to  increasing  the  voltage,  ^ n,  and  spacing  simultan¬ 


eously  to  maintain  zero  field  at  the  emitter. 
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The  Region  0.786<8<1.000 

In  this  region,  the  sequence  shown  in  Figure  2  is  Slightly  altered, 
as  Shown  in  Figure  1.  For  positive  the  sequence  is  the  same  as  in 

Figure  2.  As  shown  in  Figure  1,  the  decrease  of  T1^  through  negative 
values  leads  to  the  curve  E( emitter)  =  0,  curve  7  of  Figure  8.  Then  pro** 
deeding  along  the  "1/8"  curve  in  Figure  1  from  curve  8  in  the  first  quadrant 
the  sequence  runs  through  the  set  of  curves  shown  in  Figure  4.  Curves  9, 

10  and  11  in  Figure  4  are  not  found  in  the  sequence  for  0<9<0.786,  and 
exist  because  the  "1/9"  curve  enters  the  T!  >^min  region  as  shown  in  Figure  1 


Figure  4.  Sequence  of  Potential  Distributions 
for  0.786<9<1.000 

For  a  given  temperature,  the  volt-ampere  curve  for  a  given  9  is  easily 

extracted  since  the  values  of  It  and  ft  .  determine  the  ion  and  electron 

'  '  -  ■  m  '  nn  '  ' 

currents  that  pass  to  the  collector.  It  is  important  to  note  that  there 
is  a  solution  and  thus  a  current  for  all  potentials.  The  load  line  will 
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therefore  always  intersect  the  diode  characteristic  to  yield  an  operating 

point.  Calculations  indicate  that  similar  steady  state  solutions  exist 

for  finite  spacings.  Hernqvist^  and  Johnson^  have  claimed  to  the  Contrary 

that  Steady  state  solutions,  which  must  also  satisfy  a  load  line ,  do  not 

exist  for  certain  values  of  I.  as  will  be  pointed  out  later  it  is  actually 

possible  to  have  two  different  operating  points  with  different  currents  at 

the  same  potential  difference.  The  explanation  for  the  oscillations  (in  — 

time)  appears  more  likely  to  be  found  in  the  instability  of  the  steady 

2. 

state  solutions  in  a  manner  discussed  qualitatively  by  Hernqvist  and 
Johnson^.  A  good  discussion  on  oscillations  is  presented  by  Auer  in 
Reference  1. 

Existence  Froofs 

It  has  been  proven  analytically  that  the  potential  solutions  exist  for 

* 

the  region  bounded  by  T1 .  =fl  andT'  =T)  (T _ ).  That  is,  the  expressions  for 

m  mm  m  m  mm 

never  become  negative.  Specific  Solutions  have  been  found  in  all  the 
other  regions.  Existence  proofs  for  the  other  regions  were  not  attempted 
since  it  appears  certain  that  the  solutions  exist,  especially  since  the 
identification  of  the  E(emitter)  *  0  lines  in  the  two  quadrants  gives  the 
results  a  certain  unity. 

SOLUTIONS  FOR  FINITE  SPACING 

The  oscillatory  solutions  for  infinite  spacing  are  equally  solutions 
for  finite  spacing  if  the  collector  is  placed  at  a  position  where  all  the 
ions  and  electrons  are  moving  away  from  the  emitter.  When  this  is  not  the 
case,  the  solution  will  have  one  of  the  forms  (or  the  negative  thereof) 

shown  in  Figure  5«  Solutions  are  obtained  by  varying  8 ,  T]  and  T1  .  in 

'  .  -  ~  ~  ~  .  m  mm 

order  to  satisfy  the  boundary  conditions.  The  first  integrations  are  all 
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Figure  5*  Non-oscillatory  Potential  Distributions 
For  Finite  Spacing 

contained  in  Equations  7  to  10.  For  finite  spacing  there  are  no  equations 

Such  as  13  and  14.  As  an  example,  Equation  12  holds  for  curve  1  in  Figure 

5,  if  is  interpreted  as  the  collector  potential.  For  given  3,  Tl  , 

and  spacing  5^,  (dVdf  )a  ^_-p  would  be  varied  until  the  second  integration 

yielded  5(1  .)  =  I  . 

min  ■  L 

Of  special  interest  are  those  solutions  which  have  vanishing  electric 
field  at  either  the  emitter  or  collector,  and  are  thus  just  space-charge 
neutralized.  These  solutions  have  been  discussed  in  detail  in  the  section 


"Emission  Requirements  for  Removal  of  Space  Charge  Barriers," 

In  Figure  6,  some  solutions  for  3  =  o.i  are  shown.  It  is  seen  that 


two  solutions  with  different  currents  are  possible  at  a  given  voltage, 
since  the  two  oscillatory  solutions  are  good  finite  solutions  if  the 
collector  is  put  at  either  intersection  points  in  Figure  6.  This  bi» 
stability  has  been  previously  pointed  out  by  Auer.^  For  the  example 
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3=0.1  the  volt-ampere  curve  appears  to  have  the  shape  shown  in  Figure  7a 
rather  than  that  in  Figure  7b  which  is  taken  from  Auer's  paper.  Figure  7a 
for  a  given  3=0.1,  both  modes  correspond  to  Oscillatory  solutions  With 
potential  maximum  next  to  the  emitter.  Auer's  two  modes  in  Figure  7b  and 
his  Figure  3a  appear  to  differ  in  that  one  mode  has  a  potential  maximum 
and  the  other  a  potential  minimum  next  to  the  emitter.  For  given  B,  such 
a  possibility  does  not  seem  to  exist  according  to  Figure  3*  Such  a  possi¬ 
bility  might  exist  for  varying  .  Such  variations  in  B  may  arise  from  Ion 
barrier  charges  which  affect  the  total  arrival  rate. 
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V.  EMISSION  REQUIREMENTS  FOR  REMOVAL  OF  SPACE  CHARGE  BARRIERS 

L.  K.  Hansen  and  C.  Warnar,  III 


INTRODUCTION 

The  potential  distributions  for  CollisiOnless ,  plane  electrode ,  con¬ 
verters  have  been  examined  to  determine  the  required  conditions  for  con¬ 
verter  operation  in  the  retarding  field,  space  charge  limited ,  or  saturated 
emission  mode.  Such  analyses  have  previously  been  given  for  vacuum  diodes. 
The  following  analysis  represents  a  generalization  to  the  case  where  ions 
as  well  as  electrons  are  generated  at  the  emitter. 

VACUUM  DiODiS 

For  a  vacuum  diode  with  a  specific  spacing  the  boundary  in  volt-ampere 
space  which  separates  retarding  field  operation  from  space  charge  limited 
Operation  is  the  locus  of  points  whose  volt-ampere  coordinates  are  associ¬ 
ated  with  a  zero  field  condition  at  the  collector.  Approximate  expressions 
for  this  boundary  Curve  have  been  obtained.1-^  Similarly  the  curve  sepa¬ 
rating  the  regions  of  space  charge  limited  operation  and  saturated  emission 
operation  defines  the  parameters  associated  with  a  zero  field  condition  at 

the  emitter.  An  approximation  of  this  relationship  is  the  well  known 

4-5 

Child-Langmuir  law.  The  three  regions  of  vacuum  diode  operation  have 

also  been  defined  by  Ferris J  in  terms  of  universal  diode  characteristics. 

Nottingham^  has  used  the  same  approach.  In  his  treatment ,  however,  the 

space  charge  limited  region  appears  as  the  "collapsed"  region  between  his 

asymptotic  master  curve  and  his  universal  limiting  curve. 

7 

Langmuir  was  able  to  give  a  universal  potential  distribution  for  a 
plane  parallel,  vacuum  diode  operating  in  the  space  charge  limited  mode. 
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8 

(His  analysis  has  been  reviewed  and  amplified  by  Foritsky.  )  Langmuir 
defined  the  reduced  variables* 


1/4 


m 


(x-xm) 

m 


a) 


...  (2) 


where  L  is  the  current  density  through  the  diode ,  x^  the  position  of  the 

potential  minimum  and  V  the  value  Of  the  potential  minimum.  Langmuir's 

solution  for  x>x  =  0  is  the  zero  emitter  field  solution*  The  solution  for 

m 

x<x  =  0  is  the  zero  collector  field  solution*  Since  these  are  universal 
m 

curves  they  also  represent  the  boundaries,  in  ^-5  space,  between  the  regions 
of  retarding  field,  space  charge  limited  and  saturated  emission  operation. 

In  this  interpretation  I  can  be  thought  of  as  a  current  coordinate  (tMS?  ) 
with  a  coordinate  scale  which  depends  upon  spacing.  The  coordinate  T 
becomes  the  potential  difference  \  across  the  diode.  These  universal 
boundaries  separating  the  three  regions  of  diode  operation  are  shown  in 
Figure  1.  (Data  for  these  curves  can  be  obtained  from  Kleynen.^)  A  sample 
volt-ampere  characteristic  is  shown  by  the  dashed  line. 

The  magnitude  of  the  space  charge  barrier  is  also  obtained  easily  from 
Langmuir ' s  solution  of  the  general  potential  distribution.  A  plot  in  Tl-f 
space  of  a  curve  representing  the  constant  barrier  height  is  obtained 
simply  by  the  appropriate  translation  of  Langmuir's  solution  for  the  poten¬ 
tial  distribution.  In  the  coordinate  system  of  Figure  1,  the  translated 
curve 

|(T|L5\  =  Const*)  =  *UT'a}  *  ?(Tl  +  V  •••  (3) 

*  These  reduced  variables  differ  from  those  used  in  the  preceding  section, 
but  are  closely  related  to  the  former  variables  because  of  symmetry 
between  ions  and  electrons  in  the  space  charge  equations.  In  the  pre¬ 
ceding  paper,  |  is  measured  in  terms  of  the  emitted  current  and  T,  is  zero 
at  the  emitter. 


AI-7979 

^5 


is  the  locus  of  parameters  D  and  I  for  which  there  is  a  constant  potential 

Jj 

difference,  1^,  between  the  potential  minimum  and  the  collector.  For 

is  the  magnitude  of  the  space  charge  barrier  as  shown  in  Figure  2. 
The  translated  curve 


f(\;\=  const.)  =  mB)  ♦  i(VV  •••  C4) 

represents  the  locus  of  parameters  for  which  the  potential  barrier 
between  the  potential  minimum  and  the  emitter  is  constant.  For  Tlj>0,T!g  is 
the  magnitude  of  the  space  charge  barrier  as  well  as  the  total  potential 
barrier.  For  1L<0,  the  total  potential  barrier  is  the  sum  of  the  diode 
potential  and  the  magnitude  of  the  space  charge  barrier  These  poten 
tial  notations  are  illustrated  in  Figure  2.  The  translated  curves  are 
Shown  in  Figure  3» 

GENERALIZATION  TO  INCLUDE  ION  EMISSION 


The  equations  necessary  to  generalize  Figure  1  to  include  plane  par? 
all el  diodes  which  emit  ions  as  well  as  electrons  at  the  emitter  have 
already  been  examined  by  a  number  of  authors.^0  These  analyses  are  based 
upon  the  assumption  that  there  are  no  collisions  in  the  diode  interelec- 
trode  space.  Poisson’s  equations,  using  this  assumption,  become 

?  =  9 e‘(7,+V  {l-erf-/(H+TlL)}  -  {l+erf  >/%} 

\>0  ...  (5) 

T! 

*  dfj  =  le"71  {l+erf^L-T|}  -  |  {l-erf  1)L>0  ...  (6) 

d|3  2 

for  the  retarding  field  and  saturated  emission  cases  respectively,  where 

is  the  potential  difference  across  the  diode  and  9*  is  the  ratio  of  the 

*  The  quantity  9  given  by  Equation  7  is  the  reciprocal  of  the  9  found  in 
the  preceding  paper.  The  above  notation  is  becoming  more  conventional 
and  is  adopted  here  because  the  present  paper  is  not  a  continuation  of 
previous  work. 
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densities  at  the  emitter  surface  of  ions  and  electrons  coming  from  the 

emitter  surface. 


e 


j+ 

j- 


(7) 


The  quantities  m+  and  m-  are  the  ion  and  electron  masses  and  J+  and  J-  are 
the  current  densities  of  the  emitted  ions  and  electrons  respectively,  in 
these  equations  and  in  the  results  which  follow  the  definition  of  fl  is 
equivalent  to  that  given  by  Langmuir ,  Equation  2.  Likewise ,  §  is  defined 
in  terms  of  the  emitted  electron  current 


For  convenience  in  applying  the  results  of  this  discussion  Equation  8  is 
plotted  in  Figure  4.  The  boundary  conditions  for  the  desired  solutions  are: 


d*n 

df 

'  =  0 
^L 

\<o 

retarding  field 
Case 

(9) 

djl 

df 

=  0 

S=o 

V 

saturated  emission 
case 

(10) 

where  t(f  )  =  tl . 

L  h 

Only  monotonic  solutions  have  been  studied  and  time  and  stability  have 


not  been  considered.  The  solution  of  Equations  5  and  6,  reduced  to  a  plot 
of  vs.  for  various  values  of  3  is  shown  in  Figures  5a  and  5b,  The 
dependence  of  f ^  on  current  is  on  the  actual  emitted  current  rather  than 
the  net  diode  current,  § ,  used  in  Figure  1*  The  two  currents  are  related 
by 


h  s  1 


V° 


...  (id 


l^.leV*  \<0 

The  9=0  curves  correspond  to  the  curves  of  Figure  1.  The  shaded  region  in 


Figures  5a  and  5b  is  the  only  region  where  9>i  is  needed  to  remove  the 
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Figure  5s»  §  Emission  Requirements  for 

Removal  of  Space  Charge  Barriers 
($L  =  0  to  400) 
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Figure  5b.  9  Emission  Requirements  for  Removal 


space  charge  barrier*  At  TL~  -0.2  and  for  large  |t  the  value  rises  to  a 

la  L 

maximum  of  3=1.27.  For  a  specific  3  value  in  Figures  5a  and  5h*  "3-1  space 
is  again  divided  into  three  regions  corresponding  to  the  three  modes  of 
diode  operation  as  in  Figure  l»  It  can  be  seen,  however,  that  as  3  in¬ 
creased  from  zero  the  region  of  space  charge  limitation  diminishes  and 
eventually  disappears. 

Some  potential  distributions  corresponding  to  the  parameters  specified 


in  Figures  5a  and  5h  are  shown  in  Figure  6.  =  ±6)  It  will  be  seen 

that  for  T  <o  as  3  increases  from  zero,  eventually  a  value  is  reached  where 

not  only  the  first  derivative  but  also  the  second  is  zero  at  the  collector. 

Let  us  designate  the  value  of  when  this  occurs  as  h,  h=h(Ti^,&).  It  is 

evident  that  the  potential  distribution  for  §£h  is  independent  of  as  long 

as  §Tah.  Under  these  restrictions ,  variation  of  ?T  merely  lengthens  or 
L  u 

shortens  a  collector  plasma  column.  For  large  1^,  therefore,  and  3 
become  independent  of  and  related  only  to  each  other.  This  asymptotic 
relationship  is 

\  ,  .-2,L 

and  is  plotted  in  Figure  7.  In  a  similar  manner ,  one  finds  that  for  Tij>0, 
as  3  increases  from  zero,  a  value  of  3  is  reached  where  the  potential  dis¬ 
tribution  develops  two  inflections  and  then,  with  very  small  changes  in  3, 
goes  to  infinity  developing  a  plasma  column  which  is  isolated  from  both 
emitter  and  collector.  In  this  asymptotic  case  the  8  corresponding  to  a 
particular  T|^  is  given  by  the  minimum  of  the  following  function  of  T. 


0("3)— 1 

~ - * . . . . .  0<  I|<  '  f 

2{l-e”  }  -  2e”  {Q(\)  -  G(T)L-T1) }  " 


...  (14) 
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Figure  6.  Potential  Distributions  with  Zero 
Collector  or  Zero  Emitter  Field 


Figure  ?•  Asymptotic  B  Values 
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Where 


The  value  of  T|  for  which  the  minimum  occurs  is  the  potential  of  the  plasma 
in  the  asymptotic  potential  distribution.  It  is  zero  for  %=0  and  ap= 
proaches  0.34  as  Tl^becomes  very  large.  The  asymptotic  function,  8  (T-^)  from 
Equation  14,  is  also  plotted  in  Figure  7.  The  physical  significance  of 
these  asymptotic  values  is  simple  to  understand.  In  retarding  field  oper¬ 
ation  the  electron  current  reaching  the  collector  is  diminished  from  that 
.  _Ti 

emitted  by  a  factor  e  L.  The  ion  emission  requirements  for  neutralization 
at  the  collector  are  therefore  diminished  accordingly  since  all  emitted 
ions  are  accelerated  to  the  collector.  In  saturated  emission  operation  ion 
emission  requirements  for  removal  of  space  charge  barriers  at  the  emitter 
are  about  half  that  needed  for  neutral  emission.  This  is  understandable 


since  emitted  ions  are  reflected  in  the  interelectrode  space  and  return  to 
contribute  again  to  neutralization  at  the  emitter.  A  value  of  9=0.5  is  not 
obtained  since  complete  neutralization  is  not  needed  for  removal  of  the 
space  charge  barrier.  Since  most  of  the  ions  will  be  reflected  when  the 
potential ,  V,  across  the  diode  is  such  that  eVaskT  this  asymptotic  value 
will  be  essentially  reached  for  ^=1.  This  effect  can  be  seen  in  Figure  7. 

The  curves  of  Figure  3  can  be  placed  in  the  coordinate  system  of 
Figures  5a  and  5b  by  using  the  following  transformation 

?L  =  ?eV2  ...  (16) 

In  Figure  8  this  transformation  is  shown  for  the  same  coordinate  range  of 

Figure  5®.  Figure  8  shows  that  for  a  given  spacing  and  electron  emission 

(§=const. )  the  space  charge  barrier  (1  for  T|  <0  and  T  for  >0)  varies 

A  L  ij  Ij 

rapidly  with  TL  for  T|  <0  and  slowly  with  H  for  T!  >0.  In  Figure  5a  it  can 

jj  L  L  L 

be  seen  that  the  ion  requirement,  8,  as  a  function  of  TL  ,  varies  slowly  for 

Ju 
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TljX)  where  the  space  Charge  barrier  varies  slowly  with  ^  and  varies 
rapidly  for  ^<0  where  the  space  charge  barrier  varies  rapidly  with 
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VI.  UNIGNITED  MODE  OF  THERMIONIC  CONVERTERS 
L.  K.  Hansen  and  H.  Heri 


INTRODUCTION 

The  volt-ampere  characteristics  of  thermionic  converters  in  the  unig¬ 
nited  mode  are  sensitive  to  changes  in  emitter  temperature ,  cesium  pressure 
and  spacing.  The  dependence  of  the  characteristic  on  these  parameters  has 
been  measured  in  order  to  study  some  of  the  basic  processes  occurring  in 
thermionic  converters.  It  is  recognized  that  a  converter  operating  in  the 
unignited  mode  is  far  from  optimized  with  respect  to  either  power  output  or 
efficiency;  however,  in  this  "off-optimum"  mode  the  processes  in  thermionic 
converters  produce  detailed  structure  in  the  volt-ampere  characteristics . 
Therefore ,  the  detailed  structure  and  the  trends  in  this  structure  can  be 
used  to  identify  and  formalize  such  basic  converter  phenomena  as  potential 
distributions ,  ion  emission,  mean  free  paths,  patch  effects,  and  ionization 
processes.  The  advantage  of  using  the  unignited  mode  for  this  study  is 
that  volt-ampere  characteristics  are  not  complicated  by  the  additional 
mechanisms  which  produce  anomalously  high  current  densities  in  the  ignited 
mode . 

APPARATUS 

The  converter  used  to  obtain  the  volt-ampere  characteristics  which 
follow  was  a  variable  spacing,  plane  parallel  diode  with  molybdenum  elec¬ 
trodes  ,  *  as  shown  in  Figures  1  and  2.  The  area  of  the  plane  surfaces  of 

2 

the  electrodes  was  2.0  cm  .  Spacing  between  the  electrodes,  which  was 
variable  from  zero  to  0.125  in.  by  means  of  a  differential  screw  assembly, 

*  This  converter  was  constructed  in  another  research  program  carried  out 
under  contract  AT( 11-1) -GEN-8  with  the  U.S.  Atomic  Energy  Commission, 
(Reference  1). 
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OlAL  INDICATOR  SHAFT— 


SPACING  DRIVE  NUT 

SPACING  ADJUSTMENT  ROD 


Figure  1.  Variable- Spacing  Converter 
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Was  measured  by  means  of  an  attached  precision  dial  indicator.  The  elec¬ 
trode  surfaces  were  adjusted  with  a  precision  jig  to  make  them  parallel 
within  0.0005  in.  The  collector  was  maintained  at  a  constant  temperature 
by  the  combined  effect  of  an  attached  heater  coil  and  an  attached  cooling 
tube.  The  latter  was  cooled  by  forced  air*  The  temperature  of  the  col¬ 
lector  was  measured  with  a  thermocouple  mounted  within  the  collector. 

Visual  appearances  of  the  discharge  could  be  determined  through  a  sapphire — 
window  opposite  the  interelectrode  space.  The  emitter  temperature  was 
measured  with  an  optical  pyrometer,  sighting  through  the  sapphire  window 
and  into  a  black  body  cavity  at  the  side  of  the  emitter.  Several  layers  of 
molybdenum  foil  (not  shown  in  Figure  l)  were  placed  around  the  curved  sur¬ 
faces  of  the  emitter  to  reduce  radiation  losses. 

The  circuit  used  to  obtain  the  volt-ampere  characteristics  is  shown  in 
Figure  3.  All  characteristics  were  taken  with  an  X-Y  recorder.  Sweep 
speeds  were  about  one  inch  per  second.  The  sweep  was  controlled  manually 

!y  an  adjustment  in  the  transistorised ,  current-controlled  load. 

OLT-AMPERE  CHARACTERISTICS 

The  volt-ampere  characteristics  of  the  unignited  mode  are  shown  in 
Figures  k  to  9.  These  figures  contain  a  large  number  of  volt-ampere 
characteristics,  Such  an  extensive  coverage  of  the  parameters  concerned 
is  very  useful  for  identifying  and  studying  some  of  the  basic  effects 
occurring  in  thermionic  converters.  For  this  reason  they  are  included  in 
this  report.  An  attempt  was  made  to  develop  a  consistency  in  temperatures 
and  scale  factors  so  that  the  trends  in  the  data  could  be  easily  determined. 
For  various  experimental  reasons  this  did  not  always  occur;  however  the 
necessary  information  is  listed  with  the  individual  characteristics.  The 
essential  trends  are  relatively  easy  to  follow.  The  origins  for  these 
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Figure  4a.  Unignited  Mode  for 
Spacing  L  =  5  mils 
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Y  *  0.1  AMPS/DIV. 
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LEGEND: 

T(E)  =  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(C)  =  Collector  Temperature 


3H 

t  »  3  MILS 
X  -  2  VOLTS/OIV. 
Y  -  0.2  AMPS/OIV. 
T(E)  =  1087'K 
T(*)  «  J23”K 
T(C)  ■»  828“K 


I'lll 


4H 

L  =  5  MILS 
X  *  2  VOLTS 'DIV. 
Y  =  0.1  AMPS/DIV. 
T(E)  =  1087*K 
T(R)  =  498®K 
T(C)  =  823°K 


J _ I _ L 


Unignited  Mode  for 
Spacing  L  *  5  mils 


LEGEND: 

T(E)  —  Emitter  Temperature 

T(R)  =  Cesium  JReservoir  Temperature 

T(C;  =  Collector  Temperature 
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Figure  j?a*  Unignited  Mode  for 
Spacing  L  =  10  mils 


LEGEND : 

T(E)  =  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(C)  =  Collector  Temperature 
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Figure  5b.  Unignited  Mode  for 
Spacing  L  =  10  mils 
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L  “  10  MILS 
X  m  5  VOLTS/DIV. 
Y  =  0.1  AMPS/DIV. 
1(E)  =  1087°K 
T(t)  =  473*  K 
T(€)  ■  81l*K 
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L  -  10  MILS 
X  =  5  VOLTS/DIV. 
V  .«  0.1  AMPS/DIV. 
T(I)  #  10i7#K 
T(t)  -  44S'K 
TiC)  m  S23*K 
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1-10  MILS 
X  -  5  VOLTS/DIV. 
T  s  0.1  AMPS/DIV 
TIE)  -  10f7*K 
T(«)  9  423*K 
T(C)  -  H0‘K 


LEGEND: 

T(E)  =  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(G)  s  Collector  Temperature 
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Figure  6a.  Unignited  Mode  for 
Spacing  L  m  15  mils 


LEGEND: 

T(E)  =  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(C)  =  Collector  Temperature 
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Figure  6b.  Unignited  Mode  fop 
Spacing  L  s  15  mils 
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l  m  15  MILS 
X  =  2  VOLTS/DIV. 
Y  =  0.1  AMPS/DIV. 
T(E)  m  30«7*K 
T(«)  «  44t'K 
TIC)  -  121  "X 


LEGEND: 

He3  a  Emitter  Temperature 

T(fi)  s  Cesium  Reservoir  Temperature 

T(C)  s  Collector  Temperature 
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Figure  7a. 


Unignited  Mode  for 
Spacing  L  *  2$  mils 
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LlGlNS: 

T^E)  =  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(C)  s  Collector  Temperature 


L  P  25  Mils 
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Figure 


LEGEND: 
T(E)  = 
T(R)  s 
T(C)  = 


7b i  Unignited  Mode  for 
Spacing  L  =  25  mils 


Emitter  Temperature 
Cesium  Reservoir  Temperature 
Collector  Temperature 
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Figure  8a.  Unignited  Mode  for 
Spacing  L  =  50  mils 


LEGEND: 

T(E)  =  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(C)  =  Collector  Temperature 


AI-7979 

70 


Figure  8b.  Unignited  Mode  for 
Spacing  L  =  50  mils 


LEGEND: 

t(e)  =  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(C)  =  Collector  Temperature 
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Figure  9a 


Unignited  Mode  for 
Spacing  L  =  92  mils 


LEGEND 
T(E)  = 
T(R)  = 
T(C)  = 
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Emitter  Temperature 
Cesium  Reservoir  Temperature 
Collector  Temperature 


I  «  92  MILS 
X  =  2  VOLTS/DI V. 
Y  m  0.1  AMPS/DIV. 
T(C)  =  168 
T(*>  «  423*K 
T(C)  *  623*  K 


Fibure  9b.  Unignited  Mode  for 
Spacing  L  s  92  mils 
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T(e)  s  Emitter  Temperature 

T(R)  =  Cesium  Reservoir  Temperature 

T(C)  =  Collector  Temperature 
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Figure  10.  A  Typical  Volt-Ampere 
Characteristic 
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characteristics  were  obtained  by  temporarily  shorting  the  inputs  of  the 
X-Y  recorder. 

VISUAL  AFPEaRiiNCES 

The  visual  appearance  of  the  discharge  is  closely  correlated  to  the 
position  of  the  operating  point  on  the  volt-ampere  characteristic.  A 
typical  volt-ampere  curve  is  shown  in  Figure  10.  With  a  diode  operating 
in  region  A  of  the  characteristic ,  either  when  there  is  an  apparent 
saturation  current  as  shown  in  Figure  10  or  no  current  at  all,  there  appears 
to  be  no  radiation  coming  from  the  interelectrode  space.  In  region  B  a 
sheet  of  glow  appears  on  the  surface  of  the  collector,  increasing  the 
brightness  as  the  current  increases.  When  a  stable  region  C  exists,  a 
ball  of  glow  appears  at  the  edge  of  the  interelectrode  space,  protruding 
out  into  the  surrounding  volume.  At  region  D  a  small  ball  of  glow  appears, 
usually  at  the  center  of  the  electrodes,  as  the  current  is  increased  the 
ball  widens  to  fill  the  interelectrode  space.  For  large  spacings  this  ball 
appears  Somewhat  as  a  hemispherical  shape  resting  on  the  collector.  Some¬ 
times  the  bail  will  have  bright  columns  in  it  about  1  mm  in  diameter, 
running  perpendicular  to  the  electrodes.  As  current  is  increased  the 
columns  have  been  noticed  to  increase  in  number.  The  electrodes  are 
polycrystalline ,  therefore,  this  is  presumably  a  direct  observation  of  a 
patch  effect, 

INTERPRETATION 

A  discussion  of  these  volt-ampere  characteristics  for  the  unignited 
mode  is  presented  in  Section  VII, 
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VII.  INTERPRETATION  OF  VOLT-AMPERE  CHARACTERISTICS 
C.  Warner,  III  and  L.  K.  Hansen 


INTRODUCTION 

Only  a  small  portion  of  converter  volt-ampere  characteristics  can  fee 
used  for  efficient  energy  conversion.  The  entire  characteristic  is  useful, 
however,  as  a  tool  for  analyzing  the  basic  processes  occurring  in  thermionic 
converters.  An  understanding  of  these  processes  is  needed  at  the  present 
time  to  assist  in  evaluating  the  efficiency  limitations  in  existing  con¬ 
verters  and  to  guide  further  developments  in  thermionic  energy  conversion. 

In  the  discussion  which  follows  several  interesting  features  of  con¬ 
verter  operation  are  examined,  particularly  the  features  of  the  unignited 
mode.  Converters  operating  in  this  mode  pass  currents  which  may  be  either 
temperature  limited,  space  charge  limited,  collector  voltage  limited,  or 
diffusion  limited.  In  addition,  ions  may  fee  generated  thermally  at  the 
emitter  surface,  by  collisions  In  the  collector  sheath,  or  by  collisions  in 
the  interelectrode  space.  All  of  these  conditions  have  their  particular 
effect  on  the  volt-ampere  characteristics.  An  analysis  of  the  character¬ 
istics  therefore  leads  to  an  insight  into  these  basic  phenomena  as  they 
occur  in  thermionic  converters.  Analysis  of  the  unignited  mode  leads  also 
to  understanding  of  the  ignited  mode  because  in  the  unignited  mode  phenomena 
associated  with  ignition  can  be  studied  at  their  inception. 

Long  before  thermionic  converters  became  of  interest  for  energy  con¬ 
version  this  same  type  of  discharge  was  studied  extensively  and  was  referred 
to  as  the  low  voltage  arc  or  the  externally  heated  hot  cathode  arc.  Some 
of  this  work  is  related  to  the  following  discussion  of  converter  character¬ 
istics  and  therefore  is  reviewed,  particularly  as  it  related  to  the  unig¬ 
nited  mode  of  thermionic  converters  and  the  transition  to  the  ignited  mode. 
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REVIEW  OF  RELATED  PHENOMENA  IN  GASEOUS  DISCHARGES 

Discharges  in  general  have  what  one  might  refer  to  as  ignited  and  un¬ 
ignited  modes.  Ignition  occurs  whenever  the  increasing  current  of  a  dis¬ 
charge  introduces  a  new  phenomena  which  in  turn  increases  the  efficiency  of 
the  discharge.  Since  the  increase  in  efficiency  results  in  further  increase 
in  the  current ,  a  positive  feedback  develops  which  drives  the  discharge  into 
an  ignited  mode.  The  effect  of  the  increase  in  efficiency  is  not  only  ib- 
creased  current ,  but  also  a  reduction  in  the  applied  voltage  required  to 
sustain  the  discharge. 

In  the  Townsend  discharge ,  for  example,  positive  feedback  occurs  when 
the  ionization  level  is  sufficient  to  develop  a  sheath  with  dimensions  of 
the  order  of  the  interelectrode  space,  for  then  ion  production  becomes  more 
efficient  and  there  is  a  resulting  increase  in  secondary  electron  production. 
The  discharge  therefore  ignites  becoming  a  glow  discharge.  On  the  other 
hand,  with  the  glow  discharge ,  when  ion  bombardment  of  the  cathode  heats  the 
Cathode  to  a  point  where  thermionic  emission  supplements  secondary  electron 
emission  or  when  sheaths  develop  which  introduce  field  emission,  another 
feedback  develops  which  drives  the  discharge  into  arc  operation.  The  Volt® 
ampere  characteristic  associated  with  these  ignitions  are  shown  schemat¬ 
ically  in  Figure  1. 

Externally  heated  hot  cathode  discharges  differ  in  several  respects 
from  the  above  discharges.  In  the  ignited  or  arc  mode  a  very  striking 
difference  is  that  its  maintenance  potential  can  be  lower  than  the  ioni¬ 
zation  potential  of  the  gas,  a  feature  which  attracted  interest  almost  50 
years  ago.  Today,  of  course ,  thermionic  energy  conversion  exists  because 
the  maintenance  potential  may  be  negative .  An  explanation  for  the  ioni¬ 
zation  in  this  low  voltage  mode  was  first  demonstrated  by  Compton  and 
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GLOW  MODE 


Figure  1.  Volt-Ampere  Characteristics 
for  a  "Cold”  Electrode  Discharge 

Eckart"1"  when  they  showed  the  existence  of  a  virtual  anode  or  high  potential 

region  in  the  interelectrode  space.  Electrons  would  tend  to  neutralize  such 

2 

a  region;  however,  Druyvesteyn  proposed  an  electron-electron  collision 
mechanism  by  which  neutralization  of  this  virtual  anode  could  be  prevented. 
Johnson^  has  given  an  extensive  study  of  these  effects.  More  recently  this 
ignited  mode  has  been  studied  by  a  great  number  of  people  for  the  case  of 
cesium  in  thermionic  converters. 

For  the  discussion  which  follows ,  previous  work  with  the  unignited 
mode  of  externally  heated  hot  cathode  discharges  is  of  special  interest. 

The  unignited  mode  of  these  discharges  can  be  expected  to  differ  greatly 
from  the  unignited  modes  of  discharges  without  externally  heated  cathodes 
because  thermionic  electrons  are  available  to  the  discharge ,  and  therefore 
the  effects  of  secondary  electrons  are  suppressed.  Many  of  the  features  of 
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the  unignited  mode  presented  in  Section  VI  of  this  report  have  been  observed 

in  earlier  work  with  gases  other  than  cesium.  The  anode  glow  mode ,  for 

example ,  has  been  known  for  many  years.  An  extensive  study  of  this  mode  has 

been  made  by  Kalter ,  Johnson  and  Webster  for  noble  gases  in  cylindrical 
_  7 

diodes.  Daly  and  Emeleus  used  mercury  vapor  in  cylindrical  diodes  and 

obtained  similar  results.  A  planar  system  with  argon  used  by  Cairns  and 

8 

McCullagh  also  showed  an  anode  glow  mode.  The  intermediate  or  first 

ignited  mode  shown  by  region  C  of  Figure  10  in  Section  VI  (p  )  has  also 

.  _  .  _  _  .  .  _  q 

been  observed  previously.  Both  Daly  and  Emeleus  and  Pengelly  and  Wright 

have  examined  this  feature.  The  apparent  saturation  of  the  current  in  the 

unignited  mode  was  reported  by  Cairns  and  McCullagh  and ,  much  earlier,  by 

Duf fendack^0  using  hydrogen.  Some  of  the  results  of  Section  VI  are  at 

variance  with  the  works  referred  to  above.  For  example,  Cairns  and 

McCullagh  report  that  the  anode  glow  mode  occurs  when  the  diode  current  is 

increasing  to  the  apparent  saturation  level,  but,  as  shown  in  Figure  10  of 

Section  VI,  the  present  work  indicated  that  the  anode  glow  mode  occurred 

(when  there  was  also  an  apparent  saturation  current)  when  the  current  in- 

7 

creased  abruptly  above  the  apparent  saturation  level.  Daly  and  Emeleus 
also  ascribe  the  cause  of  their  first  ignited  mode  to  the  voltage  drop  along 
their  cathode.  They  do ,  however ,  mention  that  similar  results  were  obtained 
by  Pengelly  and  Wright  for  nearly  equipotential  cathodes.  In  the  present 
work,  no  voltage  drops  existed  across  the  cathode ;  therefore ,  the  character¬ 
istics  of  the  first  ignited  mode  do  not  appear  to  be  related  to  cathode 
voltage  drops.  This  mode,  in  which  the  glow  positions  itself  at  the  edges 
of  the  interelectrode  space,  evidently  occurs  because  the  increased  electron 
path  at  the  electrode  edges  allows  the  discharge  to  ignite  at  a  lower 
current  than  would  be  required  for  ignition  in  the  interelectrode  space. 
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The  characteristics  of  Section  VI  show  that  this  mode  occurs  Only  for  low 

cesium  pressures  where  the  extra  path  length  is  essential.  Such  an  effect 

is  common  in  other  types  of  discharges. 

AS  yet  the  unignited  mode  in  thermionic  converters  has  not  been  studied 
_  li-12 

extensively.  Houston  and  Gibbons  however,  have  measured  some  trends 
in  this  mode  for  the  case  where  there  are  many  mean  free  paths  in  the  inter* 
electrode  separation  and  have  concluded  that  the  apparent  saturation  current 
is  space  charge  limited. 

INTERPRETATION  OF  THE  UNIGNITED  MODE 


Before  the  unignited  mode  volt-ampere  characteristics  cont ained  in 
Section  VI can  be  interpreted,  several  important  parameters  must  be  calcu¬ 
lated  for  the  ranges  of  emitter  temperature  (T^)  and  cesium  reservoir  tem¬ 
perature  (T^  )  covered  by  the  characteristics.  First ,  it  is  important  to 

know  the  work  function  (0-,)  and  saturation  electron  current  ( J  )  for  the 

b  s 

molybdenum  emitter  as  a  function  of  T£  and  TR  .  This  information  has  been 

12  13 

reported  by  Houston  and  Aamodt  and  is  plotted  in  Figure  2  for  the  tem¬ 
perature  ranges  of  interest.  In  the  preparation  of  this  figure,  some  ccm- 

l4 

parisons  were  made  to  the  plots  prepared  by  Nottingham.  Some  curves  of 
constant  Debye  length  ,  h 


h  = 


J 


n — r~~» 

5  n  cm 
10  -o 


where 


*  4.02  x  10 


,13 


cm 


-3 


...  (1) 


(2) 


have  also  been  plotted  in  Figure  2.  For  convenience,  the  parameters  Tg  and 


Tg  associated  with  the  characteristics  of  Section  VI are  indicated  on  the 


figure  by  dots  adjacent  to  the  number-letter  designation  for  the  appropriate 
characteristic .  It  is  also  important  to  know  the  ratio,  at  the  emitter,  of 


AI-7979 

80 


the  ion  to  electron  densities  (0)  for  the  emitted  particles.  Lines  of  con¬ 
stant  B  have  been  plotted  in  Figure  3.  These  lines  are  calculated  from  the 
work  function  information  in  Figure  2,  Figure  3  also  contains  a  plot  of 
lines  of  constant  electron  mean  free  path*  X  ,  Nottingham's^  choice  for 
the  mean  free  path  has  been  adopted  here. 

From  Figures  2  and  3  several  limiting  cases  can  be  identified  which 
suggest  appropriate  analytical  models  for  analyzing  the  volt-ampere  charac¬ 
teristics.  Figure  3 §  for  example,  indicates  two  limiting  eases,  where  the 
mean  free  path  is  either  large  or  small  compared  to  the  interelectrode 
spacing.  In  the  latter  case,  the  surface  ionization  mode  of  the  volt- 
ampere  characteristics  are  appropriately  interpreted  with  space  charge 
analysis.  Figure  3  also  indicates  that  the  emission  may  be  either  electron 
rich  or  ion-rich  in  this  large  mean- free -path  case.  For  the  very  electron- 
rich  and  large  mean-free-path  case  the  space  charge  analysis  of  Langmuir1^ 
may  be  used  for  interpreting  the  characteristics.  If  ion  emission  is  not 
negligible,  however,  then  a  more  general  analysis  is  required,  along  the 
lines  indicated  in  Sections  IV  and  V. 

Figure  2  indicates  that  those  characteristics  associated  with  a  mean 
free  path  which  is  small  compared  to  the  interelectrode  distance  also  tend 
to  have  a  Debye  length  small  compared  to  the  interelectrode  spacing. 

Because  of  the  smallness  of  the  mean  free  path  the  analysis  of  these  char¬ 
acteristics  must  include  diffusion  and  mobility  effects.  Because  of  the 
smallness  of  the  Debye  length  a  first  analysis  may  assume  quasineutrality 
and  divide  the  interelectrode  space  into  sheath  and  plasma  regions.  Such 
a  model  is  developed  here  for  both  the  electron-rich  and  ion-rich  cases. 

The  effects  of  sheath  and  volume  ionization  can  be  identified  in  the 
characteristics  and  are  briefly  discussed  here.  An  extensive  analysis  of 
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Figure  3,  Emission  Ratio  for  Cesium-Covered1  Molybdenum 
and  Mean  Free  Path  for  Electrons  in  Cesium 


these  effects  however  will  not  he  given. 

Surface  Ionization  Ko des 
1.  Space  Charge  Mode 

A  cellisioniess  model  is  appropriate  for  interpreting  the  character¬ 
istics  of  Section  VI  which  are  associated  with  low  pressure  (Tg  <480°K)  and 
close  spacing  (L=5  mils) .  In  these  cases  the  electron  mean  free  path  is 
"large  compared  to  the  interelectrode  spacing  (Xg=2L  to  10L) .  The  collision¬ 
less  model  for  the  surface  ionization  mode  has  been  examined  in  Section  IV 
and  V.  The  analysis  of  Section  IV  provides  the  necessary  background  for 
the  calculation  Of  volt -ampere  characteristics.  Section  V  indicates  the 
voltage  interval  over  which  converter  currents  are  space  charge  limited. 

When  ion  generation  at  the  emitter  can  be  neglected  ( S«l )  the  character^ 
istics  should  approximate  those  determined  from  Langmuir • s  space  charge 
analysis. 

Since  Langmuir ' s  space  charge  analysis  assumes  a  Vacuum  diode  and 
assumes  electrons  only  to  be  generated  at  the  emitter,  the  characteristics 
most  likely  to  correspond  to  his  analysis  are  those  for  postion  7H  in 
Figures  2  and  3.  The  characteristics  for  this  position,  for  several  spac- 
ings,  is  shown  in  Figure  4.  It  is  difficult  to  compare  the  voltage  and 
spacing  dependence  obtained  from  space  charge  calculations  to  that  shown  in 
this  figure  since  the  saturation  current  is  not  determined.  Position  7G 
has  a  lower  saturation  current  and  the  results  of  Section  VI indicate  that 
the  9  value  for  this  position  is  low  enough  that  the  volt-ampere  character¬ 
istic  is  not  significantly  affected  by  ion  emission.  The  calculated  and 
experimental  characteristics  for  this  position  are  shown  in  Figure  5«  The 
theoretical  curves  were  calculated  assuming  a  saturation  current  of  62.5  ma. 
The  voltage  coordinate  for  the  experimental  curves  has  been  corrected  for  a 
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Figure  4.  Volt-Ampere  Characteristics 
for  Position  7H  of  Figure  2 

0.7-volt  contact  potential  difference.  A  possible  error  in  the  position  of 
the  curves  with  respect  to  the  current  coordinate  was  not  corrected  for. 

Figure  5  indicates  a  qualitative  agreement  between  the  experimental  and  cal® 
culated  curves.  It  is  interesting  to  note  that  the  experimental  curves 
develop  a  saturation  for  both  the  close  and  the  wide  spacing  cases.  The 
former  is  evidently  the  true  saturation  current;  the  latter  "apparent"  satu¬ 
ration  could  be  either  a  saturation  induced  by  diffusion  and  mobility 
effects  such  as  discussed  in  the  next  sub-topic  (Plasma  Mode),  or  a  true 
saturation  current  from  high  work  function  patches.  Since  the  mean  free 
path  is  so  large  in  this  case  (Ae**L)  the  latter  explanation  seems  more 
reasonable.  Since  the  two  currents  differ  only  by  a  factor  of  5  to  10,  the 
space  charge  barrier  that  would  be  needed  to  produce  this  reduction  in  current 
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is  of  the  order  of  the  difference  in  patch  contact  potential  that  would  pro¬ 
duce  the  relative  saturation  currents.  It  is  reasonable ,  therefore,  that 
some  patches  Could  be  space  charge  limited  while  others  were  emitting  satu¬ 
rated  currents.  According  to  the  results  of  Section  V  the  low  work 
function  patches  would  not  have  saturated  (62.5  ma)  for  the  50-rail  spacing 
ease  until  about  36  volts  had  been  applied  across  the  diode.  For  the  same 
spacing,  however,  the  high  work  function  patches  (5-10  ma  saturation)  would 
saturate  near  zero  applied  voltage.  This  conforms  to  the  experimental 
results  shown  in  Figure  5* 

A  result  of  Section  V  which  should  be  compared  to  the  characteristics 
Of  Section  VI  is  that  for  9  values  equal  to  or  greater  than  1.2,  no  space 
charge  limitations  can  exist  in  a  diode  operating  in  the  surface  ionization 
mode.  This  effect  can  be  observed  qualitatively  in  the  characteristics  for 
the  4-B  to  4-K,  the  5- A  to  5-H,  the  6-A  to  6-H  and  the  7-B  to  7-H  sequences, 
for  the  5-  and  10-mil  spacings.  In  the  characteristics  with  9<1.2  there  is 
little  variation  with  emitter  temperature  of  the  voltage  needed  to  produce 
saturation  current.  With  R>1 . 2 ,  however,  the  voltage  needed  to  produce 
saturation  varies  rapidly  with  emitter  temperature,  This  variation  occurs 
because  the  saturation  current  changes,  therefore  and 
exact  value  of  3  for  which  the  transition  takes  place  cannot  be  determined 
because  characteristics  were  not  taken  for  smaller  intervals  of  emitter  tem¬ 
perature  and  because  of  the  uncertainties  in  B  as  a  function  of  cesium 
reservoir  and  emitter  temperatures. 

2.  Plasma  Mode 

This  discussion  will  be  concerned  with  the  high-pressure  diode  which 
has  a  plasma  region  as  well  as  emitter  and  collector  sheaths  and  which  has 
surface  ionization  as  the  only  source  of  ions.  The  collisionless,  long 


vary.  The 
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mean- f fee-path  ease,  has  been  considered  in  Sections  IV  and  V. 
a.  Plasma  Equations 

The  basic  equations  for  the  plasma  region  are  assumed  to  be 


F 

e 


Vn 


(3) 


r 

JE 

D 

P 


Vn  + 


...  (4) 


The  notation  is  conventional  and  given  in  the  glossary.  These  equations 
assume  approximate  charge  neutrality.  They  also  assume  that  the  electron 
temperature  is  that  of  the  emitter;  and  that  the  Einstein  relation, 

U  =  eD/kT ,  is  valid.  The  electron  distribution  function  will  not  be  char¬ 
acterized  by  the  emitter  temperature  if  there  is  a  large  injection  voltage. 

In  this  case,  the  electron  "temperature"  might  better  be  determined  by  an 

17 

energy  balance  equation  or  even  more  exactly  by  approximate  use  of  the 

1 8 

Boltzmann  equation.  Inelastic  collisions  will  tend  to  keep  the  electron 
"temperature"  from  increasing  strongly  with  increasing  injection  voltage. 
For  this  reason  and  because  of  the  invalidity  of  the  model  in  the  presence 
of  important  volume  ionization  effects,  the  present  model  is  considered  a 
good  first  approximation . 


Equations  3  and  4  can  be  solved  for  the  particle  density 


r  \ 


-  - 1 


n(x)  =  n(0)  -  ^ 


r  r 

-£.  +  -£ 
D  D 

•  p 


...  (5a) 


...  (5b) 


where  n(O)  is  the  particle  density  at  the  emitter  edge  of  the  plasma , 

Equations  3  and  5a  or  4  and  5a  can  be  solved  for  the  electric  field,  £, 


AI-7979 

sar 


This  electric  field  £  then  determines  the  charge  difference  e(n  -n  )  by 

P  • 

Poisson's  equation 
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...  (7) 
or 


Equation  7  must  be  used  to  check  the  assumption  of  approximate  charge 
neutrality. 

b.  Boundary  Conditions  for  ^>0,  T1q<0 

The  plasma  equation  5b  must  be  supplemented  by  boundary  conditions 
which  reflect  the  continuity  of  particle  flow.  For  a  positive  emitter 
sheath,  =  eVg/kf>0,  and  a  negative  collector  sheath,  *  eVc/kT<0,  as 
shown  in  Figure  6,  these  boundary  conditions  are 
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...  (8) 

...  (9) 
...  (10) 

...  (11) 


where  V  and  V  are  the  emitted  electron  and  ion  fluxes.  Equation  10 
®  P 

assumes  that  the  fraction  e  of  the  current  incident  on  the  emitter  sheath 


from  the  plasma,  F  = 


neE  Ve 


-  - -g  ,  penetrates  to  the  emitter.  Equation  10 


thus  assumes  that  the  electrons  are  characterized  by  the  emitter  temperature, 
c.  Solution  of  Equations  for  ^g>0, 

The  plasma  density  n  ^  is  related  to  the  density  n@E  by  the  plasma 
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Figure  6.  Potential  Distribution 
for  Ion-Rich  Emission  and  an 
Accelerating  Collector  Sheath 


Equation  5b*  Then,  the  four  equations  (8) ,  (9),  (10)  and  (11) 
four  unknowns  F  TL,  T  ,  and  n  „  for  an  assumed  value  for  I" 

p  £j  O  6 


determine  the 
Finally ,  the 


plasma  potential  drop,  At]  is  determined  by  integrating  the  electric  field 
given  by  Equation  6. 

Algebraic  manipulation  of  Equations  8  to  11  yields 


...  (12) 


...  (13) 


(14) 


...  (15) 

...  (16) 
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where 
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The  parameter  R  has  been  chosen  so  that  the  electric  field  £  in  the  plasma 
vanishes  for  R=1  according  to  Equation  6;  the  parameter  R  vanishes  for 
infinite  applied  field  T)c=  -  eo  according  to  Equation  14. 

Plots  of  Equations  12,  14,  15  and  16  are  presented  in  Figures  7  to  10. 


In  Figure  7»  Equation  12  is  plotted  for  This  model  predicts  the 

e 

saturation  currents  in  the  ion-rich  cases.  More  explicitly,  for  R=0,  the 

1 3 

and  determines  the  saturation 


abscissa  in  Figure  7  is  just  V  =  1^7- 


current  r  /v 


R=0 


Equation  12  holds  for  all  negative  and  thus  R 


varies  from  R=0  at  tl^sago  to  RsX^/X^  at  =  0. 

Equation  14  is  plotted  in  Figure  8  for  two  values  of  X^/X^.  Equation  15 
is  plotted  in  Figure  9  and  Equation  16  in  Figure  10.  The  parameter  R  was 
varied  only  between  0  and  1,  although  the  equations  hold  for  larger  values, 
d.  Boundary  Conditions  for  TT^X),  T^X) 

The  boundary  conditions  (8)  and  (9)  at  the  collector  must  be  changed 

eyc 

for  the  positive  collector  sheath  ^  =  kT^0  shown  in  Figure  11: 
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Figure  11.  Potential  Distribution 
for  Ion-Rich  Emission  and  a 
Retarding  Collector  Sheath 

e.  Solution  of  Equation  for  T)g>o ,  T^>o 

Algebraic  manipulation  of  Equations  10,  11,  12,  17  and  18  yields 


Equations  19  to  23  are  identical  with  Equations  11  to  l6  for  T!^=0.  The  sat¬ 
uration  ion  current  for  Tl^  =  co  (R»  co)  can  be  calculated  from  Equations  19 
and  20 
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f .  Cheek  oil  Approximate  Neutrality 
According  to  Equation  7 


•  •  • 


(26) 


The  factor,  eQkT/ea n ,  is  the  square  of  the  Debye  length*  h3 .  For  T^X), 
Equation  18  is  used  to  evaluate  n: 


For  Dq<0,  Equation  8  is  used  to  evaluate  n: 


The  check  on  approximate  neutrality  is  Completed  by  evaluating  r  or  I\ 
with  Equation  24  or  12,  and  introducing  the  corresponding  density  n  into 

G 

Equation  27a  or  27b. 

g.  Diffusion  and  Conduction  Currents 

The  total  currents  r  and  r  consist  of  diffusion  currents  r  and 
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These  currents  are  related  to  the 
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parameter  R  by  the  following  equations: 
p(e) 

»  |  (1-R)  ...  (29) 


for  electrons,  and 
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for  ions. 
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(30) 


Thus,  for  the  ion-rich  case  with  an  accelerating  potential  T1  <0 ,  the  R 
varies  between  R=0(11c=-  go)  and  R=1  (6=0),  The  electron  diffusion  current 
is  seen  to  contribute  from  50%  at  R=0  to  100%  of  the  current  at  R=l.  The 
ion  diffusion  current  exceeds  the  total  ion  current  and  is  opposed  by  the 
conduction  current,  and  exactly  canceled  for  Tl^s-  co. 
h.  An  Ion-Rich  Example 

An  example  of  an  ion-rich  plasma  mode  is  given  in  the  following  dis¬ 
cussion.  For  this  example  the  following  parameters  are  chosen: 

8  s  IQ,  \e/xp  «  2,  3i*/4X#  =  100 
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ms  12  to 

16  and 

19  to 

23  yield  the  results 

shown  in 

Table  I. 

TABLE 

I.  DIODE  PARAMETERS 

FOR 

e  =  10 

’  vs  ! 

=  2,  3L/4Xe  = 

=  100 

r  /u 

ex  e 

R 

V 

.  Je. 

V(V2) 

ATI  1 

20,060 

0 

1.20 

«  ob 

101 

-4.61 

CO 

0.058 

0.06 

1.20 

-2.84 

106 

-4.12 

8.16 

.  0.050 

0.19 

1.20 

-1.73 

_ 120 

-3.25 

-4;.  18 

0,040 

0.52 

1.20 

-0.90 

153 

-1.56 

3.66 

0.032 

0.91 

1.20 

-0.47 

182 

-0.25 

1.92 

0.025 

1.43 

1.20 

-0.18 

244 

0.98 

0.40 

0.020 

2.00 

1.20 

0.00 

301 

1.90 

-0.70 

0.020 

2.10 

1.20 

0.025 

311 

2.02 

-0.85 

0.016 

2.94 

1.20 

0.22 

396 

2.76 

-1.78 

0.013 

3.95 

1*20 

0.40 

497 

3.30 

-2.50 

0.010 

5.51 

1,20 

0.61 

640 

3.74 

-3.15 

0.004 

14.64 

1.20 

I.43 

1570 

4.69 

-4.92 

0.002 

38.9 

1.20 

2,33 

4010 

5.10 

-6.23 

0 

00 

1.20 

00 

CD 

5*31 

-  00 

-ue  of  Tig 

is  assumed  to 

be  1.20  although  it 

actually 

varies  slowly 

20  at  R=0  to  l.l6  at  1 

!=«.  The 

ratio  of  the  density 

neE  to  neC  18 

m  by  the 

argument  at  the  logarithms  in  Equations  16 

and  23  for 

R>  X  /X  .  The  density  is  thus  decreased  by  a  factor  ranging  between  101 
e  p 


P 

and  3Q1. 


This  factor  is  of  the  order  of  ?L/4X  for  many  mean-free-path 


cases.  In  order  to  proceed  further,  it  is  assumed  that  the  emitter  is 


molybdenum  and  its  temperature  is  1530°K. 
Then  according  to  Figures  2  and  3, 

X@  a  5  mils 


=  2.92  ev 

=  8  *  10"2  amps/cm9 
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The  resulting  volt-ampere  curve  is  shown  in  Figure  12.  An  important  feature 
in  Table  I  and  Figure  12  is  the  saturation  of  the  diode  current  well  below 
the  emitted  current. 


•4 . -1.0 .  "  0  . . . I . Z  3 

V  (volte) 


Figure  12.  Volt-Ampere  Characteristic 
for  Ion-Bich  Emission 
According  to  Plasma  Model 

For  T1q<0,  the  Debye  length  assumes  its  largest  value  for  R=2  and  according 
to  Equation  28b  and  Figure  2: 


For  Tq>0,  the  Debye  length  assumes  its  largest  value  for  R=2  and  the  same 
inequality  holds  as  for  ^<0.  It  is  difficult  to  determine  the  extent  which 
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the  above  breakdown  of  charge  neutrality  affects  the  volt-ampere  curve  of 
Figure  12. 

1.  Patch  Effect  for  an  Ion-Rich  Case 

19 

The  importance  of  patches  has  recently  been  emphasized.  The  patch 
effect  Can  be  included  in  an  approximate  manner  by  altering  the  expressions 


for  v  and  v  in  the  boundary  conditions.  For  example,  two 
e  p 

shown  in  Figure  13  and  the  plasma  potential  is  assumed -tb  X: 


the  sur¬ 


face  potentials  of  both  patches  as  indicated. 


EMITTER 


Figure  13.  Potential  Distribution  for  Ion-Rich 
Emission  at  a  Patchy  Emitter 

This  case  is  interesting  if  the  electron  current  is  predominantly  that  of 

patch  2,  in  which 

v  fv  _ 
e  e§ 

The  ion  current  is  here  supplied  to  the  plasma  by  both  types  of  patches  in 


proportion  to  their  area 


kT 
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The  equations  12  to  16  can  now  be  used  with  the  interpretation  for  v  and 

500  V  500  V  e 

v  just  derived.  It  must  be  remembered  that  8  a  ■  -■  P  -»  rr—  P1  ,  and 
P  »  1  *5 


j .  Boundary  Conditions  for  T^CO,  T|g<6 

Equations  8  and  9  are  valid  at  the  collector.  At  the  emitter,  the 
boundary  conditions  are 


v  e 
e 


neEve 


-  -|  =  r 

2  e 


v 

P 


neEVo 


III 

2 


=.  r 


k. 


T  2  I  p 

Solution  of  Equations  for  ^<0,  Pg^O 


Algebraic  manipulation  Of  Equations  8,  9,  31  and  32  with  33  yields 


r_  x  r  r 

P  e  e 


1.E  j* 


with 


'  3L  (1+R)  J  L 


1  forf^—  »1 

e 


>  4  [i .  <¥’  *■  ] 


....  (3*0 


The  quantities  Tl^f  and  ATI  are  given  by  Equations  14,  15  and  16.  The 
saturation  current,  F  /v  ,  is  given  by  evaluating  Equation  33  for  R  s  0 

c  8 

according  to  Equation  14: 


R  =  0 


...  (35) 
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For  0 =  2,  the  ion  current  F  contributes  negligibly  to  the  total 

P  P 

current  =  F  ~r  . 

e  p 

1.  Boundary  Conditions  for  Tg<0,  ^>0 

For  T|^<o  and  *n^>0 ,  the  boundary  conditions  are  given  by  Equations  17, 
18,  31  and  32. 

li  Solution  of  Equations  for  Tl^O,  T^,>0 

Algebraic  manipulation  of  Equations  17,  18,  31  and  32  with  38  yields 


F 

v 


with 


F  * 


r 


e 

v 

e 


'e 


The  quantities  TL,  F  and  ATI  are  given  by  Equations  21,  22  and  23* 

v  r 

n.  An  Electron-Rich  Example 

An  example  of  an  electron-rich  example  is  given  in  the  following 

cuasipn.  For  this  example  the  following  parameters  are  chosen: 

1  .  3L  *<L 


dis- 


B  = 


100 


=  100 


=  2 


Equations  14,  15,  16,  33  and  34  for  ^^>0  and  Equations  21,  22,  23*  36  and 
37  for  “H^O  yield  the  results  shown  in  Table  !!• 
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TABLE  II.  DIODE  PARAMETERS  FOR 


* s  ite>  lb  *  W.  VS  * 4 

e 


r  /v 

e'e 

R 

Je 

% 

At 

t=nE-At-tc 

0.00200 

0 

-2.28 

-  <» 

-4.61 

00 

0.00186 

0.06 

-2.28 

-2.85 

-4.16 

4*73 

0.00170 

0.19 

-2.28 

-1.73 

-3.25 

2.70 

0.00132 

0.52 

-2.28 

-0.90 

-1.56 

O.18 

0.00105 

0.91 

-2.28 

-0.47 

-O.25 

-1.56  • 

0.00083 

1.43 

-2.28 

-0.18 

0.98 

-3.08 

O.OOO68 

2.00 

-2.28 

0.00 

1.90 

-4.18 

O.OOO65 

2.10 

-2.28 

0*03 

2.02 

-4*33 

0.00050 

2.94 

-2.28 

0.22 

2*76 

-5.26 

0.00040 

3.95 

-2.28 

0*40 

3.3O 

-5.98 

0.00032 

5-37 

-2.28 

0.6l 

3.74 

-6.63 

0.00013 

14.65 

-2.28 

1.43 

4.69 

— 8.4o 

0.00005 

38.9 

-2.28 

2.33 

5.10 

-9.71 

0 

00 

-2*28 

00 

5.31 

-  00 

The  volt-ampere  curve  is  plotted  in  Figure  14  for  an  emitter  temperature 


Figure  14.  Volt-Ampere  Characteristic 
for  Electron-Rich  Emission 
According  to  plasma  Model 
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of  1530°K.  An  important  feature  of  Table  II  and  Figure  14  is  the  saturation 
of  the  diode  current  at  a  value  well  below  the  emitted  current.  The  magni¬ 
tude  of  this  saturation  current  is  not  determined  primarily  by  space  charge 
limitation  but  rather  the  back  scattering  of  electrons^® by  the  gas. 

The  check  on  neutrality  with  Equation  28a  yields 


n  -a 

,e  V. 


s  16 


This  breakdown  of  the  key  assumption  on  charge  neutrality  means  that  the 

results  must  be  interpreted  only  qualitatively. 

o.  Comparison  of  the  Theory  with  Some  Experimental  Saturation 
Currents 


Equation  35  for  the  saturation  current  in  electron-rich  cases  was 
checked  for  TR  »  573°K,  Tg=  1293°K  (position  IF,  plotted  in  Figure  15)  and 
for  Tr  =  573°K,  1^  =  l422°K  (that  is,  position  IE) .  In  Figure  16,  the  ex¬ 
perimental  saturation  currents  are  plotted  versus  the  inverse  spacing 
(1/L) .  When  Tfi  and  Tg(and  therefore  also  9)  are  given,  the  points  lie  on 
a  straight  line  in  accordance  with  Equation  35*  The  saturation  current 
predicted  by  Equation  33  is  actually  a  factor  of  four  too  small ,  if  B,  X  and 

v  are  taken  from  Figures  2  and  3.  The  dependence  of  the  saturation  current 

17,22 

on  spacing  is  not  found  in  theories  with  uniform  plasmas. 

For  the  very  ion-rich  case  TR  =  423°K,  Tg=  1552°K  (that  is,  position 
7D  in  Figures  2  and  3)  the  saturation  current  is  relatively  constant. 
Equation  12  predicts  that  the  saturation  current  is  the  emitted  current  Vg 
for  V  £  1C'5  in  the  many  mean-free-path  case.  For  the  case  at  hand, 

W\e  s  1  and  the  emitted  current  would  still  be  expected.  The  observed 
current  •*  50  ma  is  actually  about  seventy  times  the  predicted  emission 
current  by  Figure  2.  High  work  function  patches  may  be  the  reason  for  this 
anomaly. 
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Sheath  Ionization  Mode 


As  the  collector  voltage  is  made  increasingly  negative  in  the  diffusion- 
limited  surface  ionization  mode ,  the  collector  sheath  voltage  finally  ex¬ 
ceeds  the  cesium  ionization  potential.  Cesium  ions  produced  in  the  collec¬ 
tor  sheath  are  able  to  neutralize  more  electrons  and  the  current  increases 


above  the  plateau  of  the  saturation  current.  This  mode  has  been  studied  by 
W.M.  Webster,  E.O.  Johnson  and  L.  Malter^’^  but  not  for  the  planar,  close- 
spaced  geometry.  The  theory  already  presented  for  the  diffusion-limited 
surface  ionization  mode  is  easily  altered  to  include  sheath  ionization  by 
use  of  Webster,  Johnson  and  Malter's  ion  production  term  in  the  boundary 
conditions.  For  example,  Equation  7  would  be  modified  to 


*.CVp 


*1**1 


a  p 


/cT 


,AV  \3 


where  C  is  defined  by  the  Child-Langmuir  law 

r  *  o 

© 

and  a  by  N  s  a  p  (V-V^) 

where  N  is  the  number  of  ions  produced  by  centimeter  of  path  lengths  per 
electron  in  a  gas  at  pressure  p  and  with  ionization  potential  V.  The  ad¬ 
ditional  electron  current  due  to  electrons  produced  in  the  sheath  is 
neglected  here,  but  could  be  included. 

In  some  cases,  the  initial  departure  of  the  electron  current  from  the 
plateau  may  be  due  to  the  increase  in  thickness  of  the  collector  sheath  with 
accelerating  voltage.  The  corresponding  decrease  in  the  plasma  thickness 
results  in  less  back-scattering  and  increased  current.  An  example  of  the 
sheath  ionization  mode  is  found  in  volt-ampere  curve  4f  for  92  mils  spacing. 
The  collector  glow  is  observed  at  those  points  where  the  current  has  broken 
away  from  the  plateau  just  prior  to  ignition. 
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In  some  cases  the  ignition  contact  potential  difference  in  the  Sheath 
ionization  mode  is  less  than  the  ionization  potential  of  cesium.  A  possible 
explanation  is  that  the  diode  emitter  has  some  high  work  function  paths  in 
Some  unexplained  way.  Calculations  Indicate  that  the  electric  field  of  the 
collector  sheath  is  not  sufficiently  Strong  to  lower  the  ionization  poten¬ 
tial  appreciably  and  this  possibility  was  discarded.  More  likely,  a  two- 
step  ionization  process  is  important ,  according  to  some  simple  calculations , 
but  the  density  of  excited  states  is  difficult  to  estimate. 

A  systematic  study  of  the  volt-ampere  curves  of  Section  IV  has  not  yet 
been  completed;  but  the  following  points  may  be  noted.  In  the  sequence  IF, 
2F,  4F,  5F,  6F  and  7F  at  $0  mils  the  saturation  current  increases  from 

2  ma  to  30  ma,  whereas  the  emitter  current  decreases  from  40  amperes  to 

3  ma.  These  saturation  currents  are  interpreted  as  diffusion-limited 
according  to  the  plasma  model  for  IF  through  6F.  For  these  electron-rich 
cases  the  emitter  sheath  also  limits  the  current.  For  7F,  which  is  ion- 
rich,  the  magnitude  of  the  current  indicates  that  the  full  emitted  current 
is  being  passed.  The  factor  of  10  between  the  observed  30  ma  and  the 
expected  ma  is  probably  associated  with  a  low  work  function  patch.  The 
additional  current  due  to  sheath  ionization  first  appears  in  3F.  For  IF 
and  2F,  the  current  at  the  emitter  (the  saturation  current)  is  too  small  to 
generate  appreciable  ions.  The  additional  current  due  to  sheath  ionization 
increases  in  3 F  through  6F,  This  increase  is  probably  due  to  increased 
effectiveness  of  the  sheath-produced  ion  in  the  larger  mean  free  path  case. 
The  number  of  ions  produced  in  the  sheath  is  very  likely  decreasing  in  the 
sequence  3F  to  6F  since  the  density  of  cesium  atoms  is  decreasing  more 
rapidly  than  the  saturation  current  is  increasing.  For  7F,  the  increase  in 
current  due  to  sheath  ionization  is  again  small.  The  reason  is  very  likely 
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that  all  the  electrons  are  traversing  the  diode  and  thus  additional  neutral¬ 
ization  does  not  increase  diode  current. 

Volt-ampere  curve  4F  for  the  sequence  92,  50,  25,  15,  10  and  5  mils 
shows  the  effect  of  spacing.  The  saturation  current  varies  inversely  with 
spacing  in  accordance  with  Equation  33  Tor  all  spacings  except  at  5  mils. 

At  5  mils,  the  close  spacing  removes  the  space  charge  barrier  and  there  is 
no  saturation  of  the  unignited  current .  - 

Finally,  as  mentioned  earlier  in  ion-rich  cases  in  position,  ?D  for 
example ,  the  saturation  current  is  almost  two  orders  of  magnitude  greater 
than  the  expected  current.  This  large  factor  is  a  strong  indication  of  low 
work  function  patches  due  perhaps  to  cesium  on  an  oxide  layer. 

IGNITED  NODE 

This  study  has  emphasized  the  unignited  mode,  but  a  few  points  can  be 
made  about  ignition  and  the  ignited  mode.  For  small  mean  free  paths,  the 
ignition  voltage  depends  on  spacing.  For  example,  the  ignition  Voltage  for 
IF  increases  from  1  volt  at  10  mil  spacing  to  3*7  volts  at  50  mils.  This 
effect  can  be  attributed  to  the  smaller  number  of  ions  being  generated  since 
there  is  less  current;  and  to  the  decreasing  effectiveness  of  such  ions  as 
the  spacing  is  increased. 

At  small  mean  free  paths,  such  as  70  in  Figure  5»  the  ignition  voltage 
is  independent  of  spacing.  A  simple  explanation  is  that  the  ignition  takes 
place  outside  the  interelectrode  space. 

<it  high  temperatures,  such  as  for  position  4a,  there  is  a  continuous 
transition  to  the  high  current  mode  because  neutralizing  ions  exist  in  great 
abundance . 

The  initial  currents  have  not  been  studied  because  the  lack  of  guard 
rings  makes  the  interpretation  difficult. 
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VIII.  VAPORIZATION  AND  DEPOSITION  AT  CESIUM  COVERED  SURFACES 

L.  K.  Hansen 

INTRODUCTION 

The  effect  of  adsorbed  cesium  on  electrode  work  functions  has  been  a 
subject  of  study  for  many  years.  In  recent  years  interest  has  increased 
due  to  the  application  of  these  effects  to  thermionic  energy  conversion. 

The  use  of  adsorbed  cesium  for  controlling  emitter  and  collector  work 
functions  and  consequently  for  increasing  converter  efficiencies  is  well 
understood.  The  effect  of  adsorbed  cesium  on  the  processes  involved  in 
converter  lifetime,  however,  has  not  received  such  extensive  study.  The 
purpose  of  this  investigation  was  to  examine  the  effect  of  adsorbed  cesium 
on  the  material  transport  processes  in  thermionic  converters,  in  partic¬ 
ular,  the  vaporization  and  deposition  processes.  In  order  to  accomplish 
this,  the  experiment  was  performed  in  a  cesium  atmosphere  with  cesium 
pressures  low  enough  so  that  the  experimental  results  were  not  affected  by 
scattering  and  diffusion  in  the  space  between  the  vaporizing  and  collecting 
surfaces. 

At  one  time  it  was  thought  that  adsorbed  cesium  would  lower  the  free 

energy  of  vaporization  and  there fore  increase  the  rate  of  vaporization  of 

1  2 

a  substrate  material.  On  the  other  hand,  Richardson  et  al.  have  assumed 

that  adsorbed  cesium  would  reduce  the  rate  of  vaporization,  producing  a 
(1-0)  dependence  where  0  is  the  percentage  cesium  coverage  of  the  surface. 
Unfortunately  the  results  of  their  experiment  were  not  sensitive  to  this 
effect  and  therefore  the  hypothesis  could  not  be  tested.  Some  measure¬ 
ments  were  made,  however,  of  the  diffusion  effect  in  the  cesium  gas,  an 
effect  purposely  eliminated  in  the  present  study.  A  (1-0)  dependence  has 
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previously  been  postulated  for  adsorbed  oxygen  on  molybdenum^  in  order  to 

4  5 

explain  Observed  effects  of  an  electric  field  on  Vaporization  rates*  ’ 

KINETICS  OF  VAPORIZATION  AN©  DEPOSITION6’7 

The  process  of  vaporization  is  closely  related  to  the  process  of 

g 

crystal  growth.  The  present  concept  of  both  processes  began  with  Volmer's 

suggestion  that  these  phase  changes  proceed  stepwise  through  an  inter- 

- mediate ,  self  adsorption  state.  This  hypothesis  was  suggested  by  the 

9  io 

earlier  work  of  Kossel  and  Stranski  who  showed  that  there  is  a  variation 
of  binding  energies  among  the  surface  atoms  of  a  crystalline  solid.  Also, 
the  existence  of  the  self  adsorbed  layer  has  been  indicated  by  the  experi¬ 
ments  of  Langmuir,11  Knudsen,1^  and  Volmer.1^  Kossel' s  idealized  model  of 
a  low-index  plane  is  shown  in  Figure  1.  In  this  figure  atoms  in  the  self 
adsorbed  state  are  marked  C. 


Figure  1.  Kossel 's  Idealized  Model 
of  a  Low- Index  Plane 
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The  vaporization  process  is  thought  to  consist  first  Of  dissociation 

from  monatomic  ledges  (positions  A  to  B  and  then  to  C).  Supposedly,  the 

self  adsorbed  atoms  then  diffuse  in  areas  between  ledges  and  finally  leave 

the  surface.  This  model  for  the  kinetics  of  the  phase  Change  has  Since 

14  15 

become  highly  developed  by  Burton,  Cabrera  and  Frank  *  and  others  for 

the  process  of  crystal  growth.  Knacke ,  Stranski  and  Wolf f^  have  developed 

the  model  for  the  process  of  vaporization.  More  recently  Birth  and 
17  i3  _ .  _ 

Pound  *  have  also  advanced  the  model  for  the  vaporization  process* 
Adsorbed  Cesium  can  be  expected  to  affect  the  dissociation  from 
ledges  and  the  surface  diffusion  of  the  self  adsorbed  atoms,  but  an 
analysis  of  the  effect  of  adsorbed  cesium  on  these  vaporization  processes 

.  17 

has  not  yet  been  made.  However,  in  the  theory  of  Hirth  and  Pound ,  even 
though  the  effects  of  adsorbed  impurities  are  not  considered,  one  may 
consider  the  effect  that  adsorbed  cesium  will  have  on  the  parameters  of 
the  theory  and  thereby  obtain  an  insight  into  the  effect  of  cesium  on  the 
vaporization  rate.  Most  of  the  parameters  of  their  theory  which  could  be 
affected  by  adsorbed  cesium  occur  in  a  common  factor  3, 


Q 


°s  exp  4F| 
v  kT  I 


where  D  is  the  surface  coefficient ,  v  is  the  vibrational  frequency  of 
atoms  at  the  surface  and  4F  is  the  free  energy  change  in  the  desorption 
process.  Under  the  basic  assumptions  of  their  theory,  and  under  the 
assumption  that  ledges  reach  equilibrium  velocity  in  a  distance  small 
compared  to  crystal  dimensions ,  the  rate  of  vaporization  is  independent  of 
Q,  and  therefore  independent  of  the  effect  of  adsorbed  cesium  through  the 
above  parameters.  This  conclusion  must  not  be  considered  a  final  answer. 
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however,  for  some  of  the  assumptions  of  Hirth  and  Pound  will  have  to  be 

reconsidered  if  adsorbed  impurities  are  involved. 

The  effect  Of  adsorbed  cesium  on  the  deposition  of  atomic  beams  has 

also  not  been  reported,  information  is  available,  however,  concerning  the 

19  20 

effect  of  Surface  oxides  and  some  adsorbed  gases  other  than  cesium,  ’ 
Rapp,  Hirth  and  Pound^  suggest  four  ways  in  which  the  condensation  coef¬ 
ficient  may  be  affected  by  adsorbed  impurities :  1 )  adsorbed  impurities  may 
diminish  thermal  accommodation  of  incident  atoms  So  that  these  atoms  are 


reflected  instead  of  adsorbed;  2)  because  of  reduced  chemical  activity  of 
the  surface ,  the  desorption  energy  will  be  reduced  and  therefore  the  mean 
time  of  adsorption  will  be  less.  This  in  turn  decreases  the  probability 
that  the  adsorbed  atom  can  diffuse  to  a  position  where  it  can  be  incorpo¬ 
rated  into  the  crystal  lattice;  3)  impurities  will  preferentially  adsorb 
at  ledges.  This  may  hinder  the  incorporation  of  self  adsorbed  atoms  into 
lattice  sites  and  may  make  deposition  dependent  on  surface  nuclei  formation; 
4)  impurities  will  affect  the  energetics  of  surface  nuclei  formation.  This 
in  turn  will  change  the  point  of  critical  supersaturation. 

Although  the  above  effects  seem  certain  to  modify  the  condensation 

coefficient  Rapp,  Hirth  and  Pound*^  found  that  a  background  gas  pressure  of 
—4 

10  mm  Hg  did  not  affect  the  condensation  coefficient  of  cadmium.  How¬ 
ever,  they  used  a  very  high  supersaturation  and  a  very  low  substrate  tem¬ 
perature.  As  an  example  of  the  effect  of  adsorbed  impurities  on  deposition, 

22 

one  may  cite  the  work  of  Sears  who  has  analyzed  the  effect  of  adsorbed 
films  on  crystal  growth  kinetics  and  has  applied  this  analysis  to  chemi¬ 
sorbed  oxygen  on  potassium.  Potassium  was  expected  to  grow  by  a  screw 

dislocation  and  therefore  was  not  expected  to  display  a  critical  super- 

23 

saturation.  Hock  and  Neumann  however  found  a  critical  supersaturation. 
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Sears  was  able  to  explain  this  unexpected  result  in  terms  of  impurity  poi¬ 
soning  of  the  dislocation  ledge.  Presumably  rotation  of  the  dislocation 
could  proceed  only  when  supersaturation  was  great  enought  to  allow  nucle- 
ation  of  new  ledges  at  the  poisoned  ledges.  These  results  are  true  in 

iq 

general;  that  is,  only  On  clean,  chemically  active  surfaces  with  low 
disregistry  between  lattice  constants  can  one  expect  to  have  a  condensation 
coefficient  close  to  unity  at  elevated  temperatures  and  low  incident 
fluxes.  On  the  other  hand,  surfaces  with  oxide  layers  and/or  adsorbed 
gases  display  a  critical  super sat ur at ion  and  therefore  have  a  condensation 
coefficient  of  unity  only  at  very  low  substrate  temperatures  and  high 
inc ident  fluxes . 

EXPERIMENT 

In  order  that  the  effect  of  adsorbed  cesium  on  vaporization  and 
deposition  couid  be  studied  the  experiment  was  performed  in  a  cesium 
atmosphere.  Preventing  this  atmosphere  from  introducing  diffusion  effects 
into  the  results  of  the  experiment  placed  severe  requirements  on  the 
experimental  approach.  Cesium  pressure  was  limited  so  that  the  vaporizing 
surface  (ribbon  filament )  and  the  collecting  surface  were  within  a  mean 
free  path  of  each  other.  Because  of  this  low  pressure  the  ribbon  filament 
was  limited  in  temperature  in  order  to  retain  a  significant  cesium  coverage. 
For  the  same  reason  filament  materials  with  low  work  functions  were 
excluded  from  the  study.  It  was  also  necessary  that  the  ribbon  filament 
be  hot  enough  and  its  vapor  pressure  high  enough  so  the  vaporized  flux 
would  be  great  enough  to  do  the  experiment  within  a  reasonable  time.  In 
addition,  the  material  for  vaporization  and  an  operating  temperature  had 
to  be  chosen  so  that  cesium  would  neither  react  nor  amalgamate  with  the 
filament .  The  material  chosen  for  the  experiment  was  nickel.  Because  of 
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nickel's  relatively  high  melting  point  (1453°C),  high  work  function 
(4,6  ev) ,  Mid  relatively  high  vapor  pressure  (about  10”^  am  Hg  at  1100®C) 
as  well  as  compatibility  with  cesium,  it  seemed  to  be  the  best  available 
material  for  the  experiment.  The  temperature  for  operation  was  chosen  as 
1025eC  where  the  vaporization  rate  is  about  6  pg/min  for  the  area  of  the 
collector.  At  this  temperature  the  filament  would  become  a  third  to  a 
half  covered  before  diffusion  and  back-scatter  from  the  cesium  became  a 
problem. 

The  filament  and  collector  separation  was  chosen  to  be  ,05  inch.  To 
allow  accurate  positioning  even  with  thermal  expansion  of  the  filament,  the 
nickel  was  applied  to  a  tungsten  ribbon  by  chemical  vapor  deposition  and 
spring  loaded  in  the  apparatus  as  shown  in  Figure  2.  The  filament  was 
mounted  in  a  stainless  steel  chamber  sealed  with  copper  shear  gaskets. 
Collector  plates  were  mounted  on  a  wheel  which  could  be  rotated  by  a 
magnetic  coupling  through  the  chamber  walls.  An  ion  pump  was  mounted  to 
the  chamber  for  final  evacuation;  initial  evacuation  was  through  a  copper 
pinch-off  tabulation.  Cesium  ampules  were  inserted  into  another  pinch-off 
tubulation  and  cesium  pressure  was  controlled  by  temperature  Control  of  a 
stainless  steel  extension  from  the  chamber.  Filament  temperatures  were 
read  through  a  sapphire  window.  The  amount  of  material  deposited  on  the 
collector  plates  was  determined  by  colorimetric  techniques. 

For  a  first  series  of  runs  the  chamber  was  baked  at  about  10-^  mm  Hg 
at  350®C.  The  cesium  ampule  was  then  broken  and  the  temperature  of  the 
chamber  lowered  to  25Q9C  for  the  experiment.  Later  runs  involved  baking 
in  hydrogen  at  atmospheric  pressure,  A  residual  hydrogen  pressure  of  about 
10  ^  mm  Hg  was  left  in  the  chamber  for  these  runs.  It  is  expected  that 
cesium  provided  some  gettering  action  during  the  experiment. 
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Figure  2.  Vacuum  Chamber  Cover  Assembly 


Collector  plates  of  sapphire,  platinum,  carbon,  stainless  steel  and 
copper  were  used.  In  every  case,  using  the  various  collector  plates  and 

the  cesium  or  the  cesium  and  hydrogen  atmosphere,  the  condensation  coeffi- 

-2. 

cients  were  less  than  l. 0-2.0  x  10  and  independent  of  cesium  pressure 
(reservoir  temperature  was  varied  between  20CC  and  200eC) .  In  these  runs 
the  wheel  temperature  was  300-400°C.  The  temperature  at  the  surfaces  of 
the  collector  plates  was  not  determined.  A  series  of  copper  collectors 
were  also  used  which  were  in  the  shape  of  hollow  cavities.  These  gave  the 
same  low  value  for  the  condensation  coefficient.  The  exact  value  of  the 
condensation  coefficient  was  not  determined  in  any  of  the  above  cases  since 
the  amount  of  material  measured  was  approximately  equal  to  the  resolution 
of  the  colorimetric  analysis,  1  fig.  Two  measurements  were  made ,  however, 
which  estimate  the  actual  magnitude  of  the  condensation  coefficient.  One 
run,  with  the  hydrogen  gas  background  and  a  filament  temperature  of  11?5“C, 
deposited  5  pg  on  a  copper  collector.  This  indicates  a  condensation  coef¬ 
ficient  of  about  1,3  x  10"3. 

In  order  that  the  application  of  these  results  to  materials  and  con¬ 
ditions  found  in  thermionic  converters  could  be  tested,  a  converter  was 
examined  to  determine  the  emitter  material  transported  under  normal  operat¬ 
ing  conditions.  This  device  had  operated  for  200  hours  with  a  molybdenum 
emitter  at  l800°C,  a  cesium  reservoir  temperature  of  350  to  400° C,  a  spac¬ 
ing  of  10  mil  and  a  collector  of  stainless  steel  at  650  to  750°C,  The 
amount  of  emitter  material  deposited  was  30  fig.  If  back-scatter  and 
diffusion  effects  due  to  the  cesium  atmosphere  and  change  in  the  vapori¬ 
zation  rate  due  to  adsorbed  cesium  on  the  emitter  are  neglected,  this 
indicates  a  condensation  coefficient  of  0,5  x  10~3. 
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CONCLUSION 

The  material  transport  measurements  obtained  here  are  the  composite 
results  of  effects  occurring  at  both  the  vaporizing  and  collecting  surfaces 
Since  no  dependence  on  cesium  pressure  was  observed  with  the  nickel  meas¬ 
urements  ,  one  may  assume  that  as  far  as  material  transport  in  thermionic 
converters  is  concerned  the  effect  of  adsorbed  cesium  upon  the  vaporization 
rate  of  the  emitter  is  not  significant  compared  to  the  effect  of  the  low 
condensation  coefficient  at  the  collector.  A  low  condensation  coefficient 
can  be  expected  when  the  depositing  surface  has  a  high  temperature  and  when 
oxides  and/or  adsorbed  impurities  are  present.  Because  of  the  dominating 
effect  of  the  low  condensation  coefficient  no  information  could  be  obtained 
concerning  the  effect  of  adsorbed  cesium  on  vaporization  rates.  The  corre¬ 
lation  between  the  results  with  nickel  and  with  molybdenum  in  an  operating 
converter  seem  to  indicate  that  not  Only  is  the  effect  of  adsorbed  cesium 
on  the  vaporization  rate  not  a  controlling  factor  for  material  transport, 
but  also,  back-scatter  and  diffusion  processes  in  the  cesium  atmosphere  are 
not  of  primary  importance. 

The  transport  of  emitter  material  and  therefore  the  lifetime  of  a 
thermionic  converter  will  be  affected  primarily  by  the  low  condensation 
coefficient  for  the  emitter  material  at  the  collector  surface.  Because  of 
the  low  value  of  this  coefficient  an  equilibrium  partial  pressure  of  the 
vaporized  emitter  material  will  build  up  in  a  converter  interelectrode 
space  and  will  greatly  reduce  the  net  vaporization  of  the  emitter  material. 
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IX.  CESIUM  PURIFICATION 
R.  L.  McKisson 


INTRODUCTION 

Since  it  has  been  shown  experimentally  that  impurities  or  additives  to 
cesium  vapor  have  a  profound  effect  on  converter  performance ,  a  program  was 
instituted  to  examine  the  purity  of  the  commercially  available  cesium*  The 
goal  of  the  program  was  to  obtain  cesium  metal  of  sufficiently  high  purity 
so  that  the  effects  of  impurities  on  converter  operation  could  be  properly 
identified.  At  the  outset ,  the  cesium  analyses  quoted  by  the  commercial 
suppliers  indicated  substantial  amounts  of  impurities  in  the  metal. 

Whether  the  metal  actually  was  so  impure ,  or  whether  the  early  analyses 
were  so  uncertain,  still  remains  an  unknown.  In  any  event,  a  program  for 
obtaining  purer  cesium  was  initiated. 

CESIUM  ANALYSIS 

The  first  step  in  this  program  was  to  contact  the  commercial  producers 
of  cesium  and  to  enlist  their  aid  in  improving  their  products.  Immediate 
interest  was  evidenced,  particularly  by  the  Dow  Chemical  Company  and 
American  Potash  and  Chemical  Corporation.  Both  companies  provided  prepar¬ 
ative  and  analytical  procedures,  and  were  cooperative  in  providing  analyses 
of  their  products. 

It  was  intended  that  first  the  purest  source  of  cesium  metal  would  be 
determined  and  them,  if  required,  an  attempt  would  be  made  to  purify  the 
metal  further  at  Atomics  International.  To  this  end,  inquiries  for  carry¬ 
ing  out  analyses  of  cesium  metal  (or  of  a  cesium  halide  to  be  prepared  by 
AI  from  cesium  metal)  were  made  at  several  analytical  laboratories,  It  was 
found ,  however,  that  the  cost  of  the  requested  analyses  was  quite  high  and 
it  was  decided  to  carry  out  the  analyses  at  AI.  There  were,  however,  no 

AI-7979 

121 


pure  cesium  salts  available  commercially  from  which  analytical  Standards 
could  be  made.  Here,  again,  assistance  was  Obtained  from  American  Fotash 
and  Chemical  Corporation  in  the  form  of  about  25  grams  of  very  highly  puri 
fied  cesium  chloride  prepared  by  their  research  laboratory  staff.  This 
material  proved  to  be  excellent,  and  was  used  as  a  standard  in  the  AI 
analyses.  The  analyses  provided  by  the  suppliers  of  recently  purchased 
cesium  indicated  that  significant  improvements  in  quality  were  being  made. 

In  Table  I  the  results  are  shown  of  the  analyses  of  cesium  metal  made 
by  the  supplier.  These  analyses  show  the  continuing  improvement  in  the 
quality  of  cesium  produced  by  Dow  Chemical  Company.  The  Dow  people  have 
been  working  on  cesium  purity  for  several  years  and  their  product  is  now 
quite  pure. 

The  AI  analytical  results  for  four  cesium  samples  are  shown  in 
Table  II.  The  results  for  cesium  from  Dow  Chemical  Company  agree  with  the 
results  shown  in  Table  1.  Cesium  from  Mine  Safety  Appliance  also  Shows  a 
high  purity,  and  the  analysis  of  the  cesium  from  American  Potash  and  Chem¬ 
ical  is  quite  good  except  for  a  high  Rb  content. 
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TABLE  ] 

:  -  TYPICAL 

CESIUM  METAL 

ANALYSES 

Sunnlier 

Element 

Dow 

9-60 

Dow 

6-61 

Kawecki 

6-61 

Am.  Potash 
&  Chemical 

Dow 

9-6-61 

Cs 

99.9+% 

99.9+% 

99.9% 

99.11% 

A1 

0.0006 

0.0005 

0.0003 

B 

<0.0080 

<0.0016 

<0.0016 

Ba 

<0.0008 

<0.0008 

~  Ca  ~ 

0.0033 

0.0010 

0.01 

0.0030 

Cr 

<0.0002 

0.0005 

eu 

0.0002 

<0.0002 

0.0011 

fe 

0.0021 

0.0006 

0.0030 

K 

0.0018 

<0.0016 

0.01 

0.03 

0.0035 

Li 

<0.0016 

ND 

0.00 

<0.0016 

Mg 

0.0003 

0.0003 

Mn 

<0.0002 

0.0002 

Na 

0.0051 

<0.0016 

0.01 

0.02 

0.0073 

Ni 

0.0002 

0.0003 

0 

0.0080 

0.0036 

Pb 

<0.0008 

<0.0002 

<0.0002 

Rb 

0.029 

0.016 

0.01 

0.39 

0.0047 

Si 

0.0042 

0.0008 

5102 

0.005 

Sn 

<0.0008 

Sr 

<0.0002 

<0.0002 

Fe+Al 

0.001 

•Certified  Analysis  of  Cesium  purchased 
P.0.  N  151 IX  -  50906  H 
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TABLE  II  -  ATOMICS  INTERNATIONAL  CESIUM  ANALYSES 


Element 

American 
Potash  & 
Chemical 

Mine 

Safety 

Appliance 

Dow 

(Kawecki) 

Penn 

Rare 

Metals 

Ag 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

Al 

<0.0005 

<0.0005 

<0.0005 

0.0020 

B*  — 

-  0.0010 

—  0.0010 

—  0.0090  - 

—  0.0150 

Ba 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

Be 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

Ca 

<0.0005 

<0.0005 

0.0005 

0.0035 

Cd 

<0.0001 

0.0004 

<0.0001 

<0.0001 

Cr 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

Cu 

<0.0005 

<0.0005 

<0.0005 

<0.0005 

Fe 

<0.0005 

<0.0005 

<0.0005 

0.0020 

Mg 

<0.0005 

0.0010 

0.0007 

<0.0005 

Mn 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

Ni 

<0.0010 

<0.0010 

<0.0010 

<0.0010 

Pb 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

* 

•r|! 

0.0070 

0.0050 

0.0140 

0.0280 

Sn 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

Ti 

<0.0005 

<0.0005 

<0.0005 

<0.0005 

Sr 

<0.0001 

<0.0001 

<0.0001 

<0.0001 

Na 

0.0020 

0.0020 

<0.0010 

0.0070 

Rb 

0.0700 

<0.0100 

<0.0100 

0.0300 

K 

0.0040 

<0.0040 

<0.0040 

0.0130 

Li 

ND 

ND 

ND 

ND 

•High  Band  Si  values  are  thought  to  be  due  to  spurious 
pickup  from  glassware* 
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NOMENCLATURE 


b 


C 


D 

e 

D 

P 

D 

s 

E 


E 

o 

E 

8 

e 


F 

F' 


AF 


f 

f 


f 

"  a 


Q(T1) 


a  Langmuir  "entropy"  factor,  defined  on  page  15  ,  Section  III 
differential  ionization  constant  discussed  in  J.D.  Cobine , 
Gaseous  Conduc tors  (1st  ed. ,  New  York,  McGraw-Hill  Publishing 
Co.,  1941),  pp  79-80 

a  Langmuir  "entropy"  factor,  defined  on  page  15  ,  Section  III 

-  *5 

exclusion  area,  nr 

constant  defined  by  Equation  51*  Section  III 
electron  diffusion  coefficient 
ion  diffusion  coefficient 

self  adsorbed  atom,  surface  diffusion  coefficient 
energy  difference  between  adatom  and  adien  (Equation  29, 
Section  III) 

energy  difference  between  adatom  and  adion  at  zero  coverage 
average  energy  of  adions  due  to  dipole-dipole  interactions 
(Equation  20  in  Section  III) 
electric  intensity 
electronic  charge 

defined  by  Equation  34,  Section  VII 
defined  by  Equation  37,  Section  VII 
free  energy  of  vaporization 

ion  penetration  factor  (Equation  26,  Section  III) 
fractional  patch  area,  Section  VII 
normal  vibrational  frequency  of  adatoms 
normal  vibrational  frequency  of  adions 
defined  by  Equation  15,  Section  V 
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G ; ( T) )  first  derivation  of  Q(Ti) 


g&  atom  pair  correlation  function 

6ai  atom-ion  pair  correlation  function 

.  _ _ 

si,®i'xl  x2'  i°n  pair  correlation  function 
H  Hamiltonian 


sa  atom-atom  interaction  energy 

H  .  atom-ion  interaction  energy  - 

;H'i  Xg)  ion-ion  interaction  energy 

Hg  Hamiltonian  for  adsorbed  particles 


h  Debye  length  or  sheath  thickness 

hc  Debye  length  at  the  collector 

I  diode  current 


1 

i 


ionization  energy 

diode  current 

emitted  electron  current 


J 

J 

k 


saturation  (Richardson)  current  density 
emitted  ion  current  density 
emitted  electron  current  density 
Boltzmann  constant 


L 

M 

M 

o 

m,m 


diode  interelectrode  spacing 
dipole  moment  of  ion-image  pair 
Dipole  moment  of  ion-image  pair  at  zero 
electronic  mass 


M,m+ 

N 


N 


Nj 


mass  of  cesium  ion 
total  number  of  atoms  in  the  gas 
total  number  of  electrons  in  the  gas 
total  number  of  ions  in  the  gas 
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coverage 


N 

sa 

N  .  *' 

si,  si 

n,n(x) 

ae 

neE«  neC 


+6 


P.P 

h 

Q 

’a 

% 

R 

?id 

* 

r 

T 

Tg,  T(E) 
Tr.  TCR) 
T(C) 


total  number  of  adatoms  on  the  surface 
total  number  of  adions  on  the  surface 
density  of  ions  or  electrons  in  the  plasma 
density  of  electrons  in  the  plasma 

density  of  electrons  in  the  plasma  and  near  the  emitter/collector 
density  of  ions  in  the  plasma 

density  of  ions  emitted  at  the  emitter - - - 

■twice  the  emitted  ion  density 

density  of  electrons  emitted  at  the  emitter 

twice  the  emitted  electron  density 

pressure 

momentum 

canonical  ensemble  partition  function  in  Section  III 
vaporization  coefficient  defined  on  page  112,  Section  VIII 
classical  partition  function  for  a  one-dimensional  harmonic 
oscillator  for  the  adatoms 

classical  partition  function  for  a  one-dimensional  harmonic 
oscillator  for  the  adions 

current  ratio  coefficient  defined  on  page  91»  Section  III 

distance  between  particles  i  and  j 

hard  core  distance  of  adatoms  and  adions 

surface  area 

free  area  defined  by  Equation  19 t  Section  III 

temperature 

emitter  temperature 

cesium  reservoir  temperature 

collector  temperature 
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u(ij) 

V 


V 


v 

e 


V„ 

P 


x,x. 


m 


O' 

B 


r  ,r 

•  P 
fCe.p) 
diff 
r(e,p) 
cond 


A 


1 


I 


hard  cere  potential 

voltage 

Volume 

ionization  potential 
Voltage  across  diode 
electron  velocity 

ion  velocity  - — — — — - 

space  coordinate 

distance  from  the  emitter  at  which  the  field  first  vanishes 
defined  by  Equation  42,  Section  111 
ion  polarizability 

in  Section  IV,  density  ratio  at  the  emitter  surface  of 

electrons  to  ions  coming  from  the  emitter  surface 

in  Sections  V  and  VII,  density  ratio  at  the  emitter  surface  of 

ions  to  electrons  Coming  from  the  emitter  surface 

random  electron  flux  in  the  plasma  near  the  emitter 

net  electron  flux/ion  flux 

diffusion  current  defined  on  page  96,  Section  VII 
conduction  current  defined  on  page  96,  Section  VII 
defined  on  page  91*  Section  VII 
permitivity  of  free  space 

dimensionless  potential  defined  on  page  Jl,  Section  IV 
dimensionless  potential  across  the  diode 
dimensionless  potential  difference  between  the  potential 
minimum  and  the  collector 

dimensionless  potential  difference  between  the  potential 
minimum  and  the  collector 
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dimensionless  potential  at  which  the  electric  field  first 
Vanishes 

dimensionless  potential  of  the  collector 
dimensionless  potential  of  the  emitter  sheath 
dimensionless  potential  of  the  collector  sheath 
dimensionless  potential  where  the  electric  field  vanishes  a 

second  time  _ _ _ _ _ _ 

dimensionless  plasma  potential  drop 

defined  by  Equation  20,  section  IV 

negative  of  the  dimensionless  collector  potential 

defined  by  Equation  21  in  Section  IV 

defined  by  Equation  22  in  Section  IV 

percent  cesium  coverage  of  a  surface 

thermal  defiroglie  wave  length  defined  on  page  12,  Section  III 

electron  mean  free  path  in  cesium 

cesium  ion  mean  free  path  in  cesium 

cesium  arrival  rate 

electron  chemical  potential  in  a  gas 

electron  chemical  potential  in  a  metal 

atom  chemical  potential  in  a  gas 

ion  chemical  potential  in  a  gas 

adatom  chemical  potential 

adion  chemical  potential 

frequency  of  normal  surface  vibration  of  self  adsorbed  atoms 
atom  departure  rate 
electron  departure  rate 
ion  departure  rate 


5 


CP 

$ 


ao 


io 


'I 


erf (x) 
erfc(x) 


dimensionless  distance  (proportional  to  the  Square  root  of 
the  diode  current  and  therefore  also  a  current  coordinate) 
defined  by  Equation  2  of  Section  IV  and  Equation  1  of  Section  V 
dimensionless  distance  corresponding  to  diode  spacing  L 
(proportional  to  the  square  root  of  the  diode  current  and 
therefore  also  a  current  coordinate) 

"dimensionless  distance  corresponding  to  1) 

m 

monolayer  density  for  Cs-[w] 

electron  work  function 

energy  defined  by  Equation  18,  Section  III 

energy  defined  by  Equation  33*  Section  III 

energy  required  to  convert  adsorbed  ion  to  free  atom 

energy  of  cesium  evaporation 

collector  work  function 

emitter  work  function 

energy  of  ionic  evaporation  at  zero  coverage 
defined  on  page  22,  Section  III 
error  function 

complementary  error  function  [l-erf (x) 3 
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